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ABSTRACT 

This report treats the flow field downstream of an oblique shock 

* 

front that is perturbed by either plane sound waves or plane shear waves. 

The sound wave analysis predicts that under certain conditions there 
exists an angle of incidence which produces shock wave resonance. This 
results in a large amplification of the sound pressure and in large ampli- 
tude motion of the shock front itself. 

Experimental studies of sound wave-shock front interaction showed 
that the shock front is definitely excited by an impressed acoustic field. 
The shock front motion is periodic and at the same frequency as the 
impressed sound field. Natural pressure fluctuations in the upstream flow 
also excited the shock front and amplifications at the shock front of 22 
to 30 dB were found for frequencies from 20 to 2000 Hz. 

The shear wave analysis contains some anomalies that have not been 
reconciled. Therefore it requires critical review and, if necessary, 
revision before being used. 
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D. 


l 


E 
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6 

H 

h 

L 
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1 

1 

1 
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0 
M. 
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P, 

R 


r 

s 

T 

t 

U. 


as defined m text 

isentropic speed of sound in Region i 
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as defined in text 
constant pressure specific heat 
constant volume specific heat 
as defined in text 

amplitude of shock front displacement 

as defined in equation (41) 
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acoustic Mach number , U /a 

o 

Mach number in Region i, Vh/a.. 
pressure of gas 
gas constant 

u x /u 2 

specific entropy of gas 
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time 

unperturbed velocity component in i-direction 5 or velocity 
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W 
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X. 
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2 . 
1 


normal Mach number in Region i, If /a. 
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velocity component in i-direction 
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velocity component tangent to shock front in Region 
unperturbed velocity component in Region i 
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7 
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0 

A 

5 

X 

p 

0 
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shock front dxsplacement 

3a 

3y 

5a 

5t 

entropy-vorticity wave angle 
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p v 
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sound wave incidence angle 
wavelength 
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shock wave angle 
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refracted-sound wave angle 
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Subscripts 
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sw 

Superscripts 

' perturbed component 
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A 
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a 
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a 

function 

defined 

in 

equations 
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B 
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b 
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m 

equation 

(121) 


C 
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C' 
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D' 

function 

defined 

in 

equations 

(116), 
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d 
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eauation 
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entropy- vorticity wave 

perturbed 

shock wave 
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E 
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F 

F' 

f 
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k 

k 
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m 
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s 

P 

P 
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SPL 

T 



function defined in equation (117) 
function defined in equation (144) 
function defined in equation (118) 
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frequency 

function defined in equation (146) 

wave number of shear wave in region A 

surface wave number of pressure wave 

stream Mach number (W,/a.) 

A A 

velocity ratio across normal shock (U^/U ) 
nodal point on shock front (Figure 10) 
period 

pressure of gas 

reference pressure (0.290 x 10“® psi . ) 

coefficient defined m equation (129) 

coefficient defined in equation (152) 

sound pressure level 

absolute temperature of gas ' 

absolute temperature of gas in region A 

time 

stream velocity component in x-direction (Figure 13) 
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Mach number associated with U (U./a.) 
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stream velocity component in y-direction (Figure 13) 
nodal point speed 
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V 
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w 
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w' 
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s 

X 
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y 
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surface pressure wave speed 
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Mach number associated with W (W/a) 
stream velocity in 5 -direction (Figure 13) 
stream velocity before oblique shock (Figure 10) 
perturbation velocity component in ^-direction 
perturbation velocity component in ^-direction 
perturbation velocity component in n-direction 
perturbation velocity component m n-direction 
perturbation velocity (Figure 10) 
perturbation velocity as defined in text 
pressure wave velocity component in ^-direction 
pressure wave velocity component in n-direction 
refracted shear wave velocity component in ^-direction 
refracted shear wave velocity component in n-direction 
rectangular coordinate (Figure 13) 
rectangular coordinate (Figure 14) 
rectangular coordinate (Figure 13) 
rectangular coordinate (Figure 14) 


Greek Symbols 


6 function defined in equation (109) 

0 function defined in equations (112),' (141) 

w 

y specific heat ratio 

6 turning angle across shock 

o shear wave phase angle 

P 

6^ shear wave phase angle 



xi 


® shock 
Sp 

6x 

t 
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e 
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0 

cr 

X 

X 

sur 

4 

5, i, c A 
n 

n* 

n** 

a 

<t> 

<i> 

4 >' 

i' 


shock deflection, phase -angle 
perturbation pressure 
local shock deflection 

measure of strength of incident shear wave (equation (79)) 

rectangular coordinate (Figures 13, 14) 

shear wave incidence angle (Figure 10) 

shear wave incidence angle (Figure 14) 

critical value of 0 (equation (80)) 

wavelength of incident shear waves 

wavelength of pressure waves 

Mach angle (tan - - 1 (I/O )) 

w 

rectangular coordinate (Figures 13, 14) 
as defined in text 

coefficient defined in equation (130) 
coefficient defined m equation (153) 
angle of shock to incoming flow (Figure 10) 
refracted shear wave angle (Figure 14) 
refracted shear wave angle (Figure 14) 
pressure wave angle (Figure 14) 
pressure wave angle (Figure 14) 


Subscripts 

A 

o 

rms 


region A (upstream of shock) 
evaluated at shock (x = 0) 
root-mean- square 



CH'APTER I 


INTRODUCTION 
Statement of Intent 

The intent of this investigation is to conduct an ana- 
lytical and an experimental investigation of the flow field 
produced by the interaction of sound waves and shear waves 
with an oblique shock front. 

The analytical phase of this work deals with an exten- 
sion of the normal shock analyses to the oblique shock case; 
the experimental phase entails a pilot study of the pressure 
field associated with the sound wave shock front interaction 
phenomenon . 

Historical Background 

Although there have been analytical studies of sound 
wave and shear wave interactions with normal shocks, there 
has been no investigation of their interactions with station- 
ary oblique shock fronts. Consequently this discussion will 
be devoted to the previous research efforts involving normal 
shocks that are most germane to the proposed research. The 
classic work in this area has been done by H. S. Ribner, 

F. K. Moore and W. R. Johnson. 

Ribner’s initial study [1] was concerned with the 
— 

Numbers in brackets refer to similarly numbered references 
in the Literature Cited„ 
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convection of a pattern of vorticity through a normal shock. 

It treated the simple case of a plane sinusoidal shear wave 
of any orientation and wavelength. The interaction of such 
a wave with a normal shock front produced: (i) a refracted 

shear wave of altered wavelength and amplitude (the intensity 
of the shear wave is increased), (ii) an entropy wave which 
moves with the flow downstream of the shock, (iii) either a 
^pressure" wave or a sound wave depending on the orientation 
of the initial shear wave (the pressure wave is quite intense 
just behind the shock, but it is exponentially damped with 
distance downstream from the shock front), and (iv) a dis- 
torted shock front, A relatively weak initial shear wave was 
found to produce a surprisingly intense pressure wave or 
sound field downstream of the shock. 

Later studies by Ribner C2] and Ram and Rifoner [3] ex- 
tended the plane shear wave case to those of the passage of 
(i) both isotropic and strongly axisymmetrical initial tur- 
bulence and (ii) a columnar vortex through a normal shock 
front „ 

Moore [4], in his investigations, considered the unsteady 
oblique interaction of a normal shock wave with a plane dis- 
turbance. His analyses treated a normal shock wave traveling 
into a perfect gas which is at "rest". Solutions were given 
for three cases of interaction: (l) plane sound waves pro- 
pagating in the "rest gas" into which the shock moves, (ii) 
plane sound waves overtaking the shock front from behind, and 
(iii) stationary incompressible vorticity waves in the 


gas 
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ahead of the shock. In the first case the sound waves were 
refracted either as simple isentropic sound waves or as at- 
tenuating lsentropic pressure waves depending on the orien- 
tation of the initial sound waves. A stationary vorticity 
wave was also produced behind the shock. In the second case 
the sound waves were reflected as sound waves and a station- 
ary vorticity wave was produced. The third case was the 
same as Ribner's initial study [li. However, Moore treated 
it as a completely unsteady problem whereas Ribner transformed 
it into a steady flow problem. 

Johnson [5], in his Ph.D. dissertation, considered the 
interaction of plane and cylindrical sound waves with a 
stationary normal shock. His plane wave case was the sane as 
the first two cases of Moore’s study except that the shock 
front was stationary in Johnson’s analysis and it was moving 
in the analyses of Moore. Johnson's analysis, however, 
strengthened the mathematical foundation for the plane wave 
case in that it was not necessary to appeal to physical 
reasoning for the form of certain functions. 

Purdy [6] considered the convection of a plane sound 
wave of arbitrary orientation and wavelength through a sta- 
tionary oblique shock front. While this appears to be a 
simple extension of the studies by Moore and Johnson, it was 
found that the shock front oscillations appear to be unstable 
for certain combinations of initial Mach number, shock front 
angle and sound wave incidence angle. This investigation 
formed the basis for the research program reported herein. 
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The most recent work appearing in the literature is 
that of Lowson [71. It contains computations based on the 
methods of Ribner [2] and Moore [4] for shear wave and sound 
wave interactions with normal shock fronts. The computations 
are compared with experimental data. 

Significance and Scope 

The study of plane shear wave and plane sound wave in- 
teraction with oblique shock fronts is an abstraction of at 
least two very significant physical phenomena, namely, the 
pressure field associated with the oscillating oblique shock 
front produced by an axisymmetric step, and the pressure 
field associated with a supersonic fluid jet. 

The first phenomenon is currently a problem area for 
large launch vehicles since sizeable pressure fluctuations 
have been found to exist at interstage flares. It is quite 
possible that these fluctuations are due to the interaction 
of turbulent wakes, from the launch escape tower and vehicle 
protuberences , with the oblique shocks at the flares. It is 
also possible that in wind tunnel tests the tunnel sound field 
interaction with these shocks also contributes to the mea- 
sured pressure fluctuations <> When oblique shocks occur in 
series as they do for many current configurations the problem 
is compounded. If, for example, shear waves interact with the 
first shock then the intensified shear waves and the con- 
comitant sound waves and/or pressure waves will all/both 
interact with the next shock. It is clear that each 
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interaction could yield a fluctuating pressure field which 
would be stronger than the previous one. For the case of a 
flare the reattachment shock should oscillate with a larger 
amplitude and produce more intense pressure waves and/or 
sound waves than the separation shock. 

The second phenomenon is characterized by multiple 
shocks and therefore the reasoning given above should be 
applicable here as well. 

Although the models that are to be analyzed are much 
simplier than the actual phenomena, they should portray the 
characteristic behavior of such processes and they should 
also predict typical orders of magnitude for the pressure 
f luctuat ions . 

The research reported herein is specifically concerned 
with the following problems: 

1. The analytical prediction of the flow field down- 
stream of a stationary oblique shock front v/hen 
plane sound waves are present in the flow. The 
sound waves are to be of arbitrary orientation and 
wavelength and they have their origin upstream of 
the shock front. 

2. The analytical prediction of the flow field down- 
stream of a stationary oblique shock front for the 
convection of a plane sinusoidal shear wave of 
arbitrary orientation and wavelength through the 


shock . 
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3. The experimental investigation of the sound pressure 
field upstream and. downstream of a stationary oblique 
shock front when sound waves interact with the shock . 

In the two analytical studies the medium will be assumed 
to be an inviscid perfect gas. The experimental study will 
utilize atmospheric air as the working fluid. The specific 
objectives of this research are: 

1. To obtain mathematical expressions for the velocity 
field, the pressure field, the density field, the 
shock front displacement, and the shock front velo- 
city as functions of the upstream Mach number, the 
shock front angle, and the angle of incidence and 
wavelength of (i) plane sound waves and (ii) plane 
sinusoidal shear waves interacting with a stationary 
oblique shock front. 

2. To obtain experimental data for the acoustic pressure 
field associated with the interaction of s ound waves 
and a stationary oblique shock front, and 

3. To compare the experimental data for sound pressure 
amplification across the shock front with the values 
predicted by th.e mathematical model for sound wave 


shock front interaction.. 



CHAPTER II 


ANALYTICAL INVESTIGATIONS*' 

A. Sound Have Shock Front Interaction 

The perturbed flow field downstream of a stationary 
oblique shock front produced by the convection-propagatiori 
of plane sound waves through the shock front is to be deter- 
mined. This analysis represents an extension of the earlier 
analyses of Ribner El] and Moore [4-]; its main objective is 
to provide a simple means for computing the change in the 
acoustic properties across a specific oblique shock front. 

Math emati cal Formulation of the Problem 

To formulate a mathematically tractable problem, it is 
assumed that the acoustic pressure is sufficiently small for 
the pressure waves to travel at the isentropic speed in the 
undisturbed uniform flow. This wave speed is relative to an 
observer moving with the flow. 

Unsteady-Flow Problem . The flow field shown in Figure 
1 is unsteady as viewed by an observer on the shock, front 
since, in general, the nodal lines will either move up or down 
the shock front. The nodal speed along the shock is given by 

A ' ” 

The analysis presented for this case is a revised extended 
version of that contained in reference [6]. 



(see Figure 2) 
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V = [Mi cos a + sin 0] a, . (1) 

swyj 1 1 

There seems to be no advantage to changing the point 
of view to one in which the nodal lines are stationary since 
the resulting problem would still be unsteady. Thus the 
governing equations will be obtained for the Cartesian co- 
ordinate system shown m Figure 1. 

Governing Equations . The assumptions upon which the 
governing equations will be based are 
Assumptions 

1. Two-dimensional flow 

2. Inviscid flow 

3. Adiabatic flow 

4. No body forces 
5„ Perfect gas 

6„ Small perturbation of a uniform flow; that is, f-or 
Cartesian coordinates with corresponding velo- 

city components u^, we have 

I I 

U 1 = + U 1 s u 2 = (Figure 3); 

i i t 

p=ptp,p=p+p , T = T + T , and 

_ i 

s = s + s 
with 

Ui/Uj , u 2 /U l5 p /p, p /p, T /T 

t _ 

and s /s much less than one* 
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The governing equations are 
1. Continuity 


1 Dp_ 
p Dt 



0 


The continuity equation has a simpler form for the flow 
under consideration , since, for this flow, 

i t 

u i = U i + u i , u 2 = U2 , u 3 = 0 and U j = constant. 


Therefore , the continuity equation is 

I T 


1 Dp 
p Dt 


3u i 3 u 2 

Tx7 + 3 x 2 


0 


( 2 ) 


2 . Momentum 


3u . 3u . 

i i 

3 1 U i 3 x . 

3 


1 3 p 

p 3x . 
l 


or if the momentum equations are written in terms of the 


perturbation quantities 


III! 


[1 + 


r _ 

3u i 

3 u i 

i 

3 u i 

T 

3u i 

p 

/p ] 

at + u > 

3 x x 

+ U 1 

9 x i 

+ u 2 

3x 2 

[1 + 

1 

( _ 

3u 2 

3u 2 

i 

t 

9 u 2 

- 

9u 2 

p 

/p ] 

at + u > 

3 x j 

+ U 1 

3xi 

+ u 2 

3 x 2 


1 

P 


1p_ 

3xj 


1 3 p 
— 3 x? 

P 


Since the sound field characterizes the perturbations 
and it is assumed to be composed of simple sinusoidal waves, 
the foregoing equations will be nondimensionalized in terms 
of the wavelength A and the circular frequency w of these 
waves. The following dimensionless variables are chosen 
such that derivatives with respect to dimensionless position 
and time will be of the order of magnitude of one and the 
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dimensionless velocity components, pressure and density will 


be of the order of magnitude of one: 




A 

XI 

= M 

2m 


- 

A 

x 2 

= *2 

2m 



~ A 

t 

= t /» 

= t 



2ira 

Ul 

- Ui 

U] 

i 

~ T 


U 1 

= 

C 

o 

ts 

i 

~ t 


u 2 

= u 2 

u 02 

p 

= P P 


1 

~ t 

Uo 

p 

= P 

P "T 



a 

p 

= P P 


T 

~ * 

Uo 

P 

= P 

TP — 



a 


where 

2 2 2 
u 0 = Uoi + U 02 

_2 _ 
a = vp/P 


The momentum equations in dimensionless form are 

„r~ Uoi Uoi U 2 ^ 3u. 

[1 + MqP /p] — ■r - + 7f Uj — r— + 

u ° at u ° i axi. 

1 * 6 1/1 1 1 111 1 


U 01 U 01 



" I 

u 


3u. U 01 U o 2 ~i 

777 + ^ ~ “ 2 777 



1611 1611 1/1 1 

Represents the order of magnitude of the term or quantity in 
the equation directly above. 
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and 




T £ 

U 0 2 

3u 2 

°02 

Ux - 

*' i 
3 u 2 


[1 

+ MoP 

/pJ 

u 0 ' 

3 1 

+ Uo 

U2 

a 

3x a 

+ 

1 

6 1/1 

1 

1 

1 

1 1 

1 




f 




- » 



Uoi 

Uo 

^02 

a 

_ t 9u 2 
Uj -~- 

Sxj 

• + 

U 02 

U 02 

I u 2 

a 

3u 2 

3x 2 

P 

ie! 

3x 2 

1 

6 

1 1 


1 

S 1 

1 

1/1 

1 


where 

U 0 = amplitude of the sound particle velocit 3 f. 


U o i = U o cos 0 , 
u 02 = u 0 sin 0 > 

6 = angle between the sound wave propagation 
direction and the xj axis. 


and where the following order of magnitudes are assumed: 
ORDER (Ux/a) = 1 and ORDER (U 0 /a) = 5 << 1. 


These assumptions are consistent with supersonic flow 
linear acoustics. 

Retaining only those terms that are of the order 
magnitude of one we have, in dimensional form, 

i i 


3u l 


+ Uj 


3 u j 

3*1 


1 

P 


3p 

3x2 


and 


of 


<3a) 


and 


T 

3u 2 

+ U x 


8u 2 



9p 

3x 2 


3t 


P 


(3b) 
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3. Equation of State 

For the assumed perfect gas we h'ave 

p = pRT 

h = C T 
P 

Tds = dh - 

P 

2 

a = yRT 

The combination of these equations yields 


dp _ dp ds 



Integration between the unperturbed state and the perturbed 

A 

state gives the following approximate relation 


t i » 



p pa p 


4. Energy 


It was assumed that there is no heat transfer between 
a fluid particle and its surroundings; thus changes of state 
are adiabatic. If the shock front is excluded from the flow 
field (the shock front will serve as a boundary condition for 
the flow field), the changes of state of a fluid particle are 



and that 


dp a 1 

2 -b* 

pa pa 


dp . 


This assumes that 
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also reversible. Thus 


Ds 

Dt 


0. 


(5) 


Reduction of the Governing' Equations 


The momentum equations have already been reduced to a 
linear form by an order of magnitude analysis 0 In order to 
further simplify the governing equations the density will 
now be eliminated from the continuity equation,, 

Since the density, pressure and entropy are Eulerian 
variables , 


dp = iP- dt + 

dp = It dt + 

ds = I®. dt + 


3 P 
3x . 

l 

dx . , 

l 

3p 

3 x . 
l 

dx . and 
l 

3s 
3 X . 

dx . . 

i 


i 


Substituting the foregoing equations into the differential 
form of equation (4), dividing by dt and taking the limit 
as dt approaches zero results in the following expression 
for a fluid particle 


1 Dp 1 Dp 1 Ds 

- Dt --2 Dt ~ C Dt 

P pa P 

Ds 

However, from equation (5), = 0, and therefore' 

1 Dp 1 Dp 1 

- Dt --2 Dt ” --2 

p pa pa 


3p 

3 1 


U 1 


9 P 

3x- 


+ «2 


ip. 


3 X ; 
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We have expressed the static pressure, p, as the sum 

of a time-mean component, p, and a time-dependent component, 

r _ 

p o Since the unperturbed flow field is uniform, p is con- 

— * 3 p *0 p 

si: ant* Thus p = p + p and consequently 

T T 

_P — = — and JlE — - ~ P — , Introducing dimensionless var- 

3xi 3xi 9x2 3x2 

iables, as previously defined, the substantial derivative of 

the density becomes 

* - U n 

_1 Dp _ 2t a 0 

p * a | p a~ L 3t a 3xx 


r 3p* 


™'t 


U 1 a 
M °P 

* U 1 -f- 


1 / 1*1 1 11 
U 02 

a 9x ^ a 3x2 j 

5 1 1 6 11 


U01 3p’ + 

uj — — + 


~ > 3p * 

u 2 


If only those terms in the bracket that are of the order of 
magnitude of one are retained, then, in dimensional form, 
we have 


1 Dp 1 

- Dt 
P p a L 


+ ui 

3 1 1 3xx 


Substituting the foregoing equation and equation (3a) 
into equation (2) we obtain 
Ui 3u j 


1 3p 

--2 9 1 
pa 


-2 3t 
a 


2-i 

Ul 


- 2 

a 


9 u j 3 U 2 
+ 


3x 


1 


9X2 


( 6 ) 


"X 

This term appears in dimensional form. 
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.Finally, the governing equations have been reduced to 
a set of partial differential equations that express the' 
acoustic pressure field in terms of the velocity; that is, 


equations (3a), (3b) and (6), rewritten as 

* 


3p _ 


with 


3xi 


a P 

9x2 

3jL 

at 


dp 


= - P 


= - p 


9 u j 9 u j 

+ Ul 


at 


9 x 


1 J 


! t -» 

9 U 2 3 U 2 

_ + 0l _ 


at 


■i J 


t 


2 

i 

2 BU2 

3u i 

_ 2 

r u ii 

3uj 

Ul at 

“ s. 

1 ” 

L a J 

9xi 

a 9 x 2 


* ft' « ♦ ^ 


. dxi + dx 2 • 

9xi 1 9x 2 z 


(7a) 


(7b) 


(7c) 


( 7d ) 


Solution of the Governing Equations 


The solutions to equations (7) in the regions upstream 
and downstream of the shock front are linked to one another 
through the oblique shock front relations . These boundary 
conditions will be presented at this point in the analysis. 


Boundary Conditions « The unperturbed and perturbed 
oblique shock system is shown in Figure 4. From Figure. 4 
it follows that 


Vi = 

W i cos 

<7 , 



( 8a ) 

Ui = 

W i sin 

O, 



(8b) 

V *P 

= Vj cos 

a t Ui 

y 

» 

sin a + uii cos 

y 

(a-a ) + 





< 

u 2 i sin 

(0 - a y > 

( 8c ) 
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and 


Ujp = Uj cos - Vj sin sin (c-a^) - 


uoj cos (cr-a> ) - a. , 


C 8d ) 


Since only small perturbations are being considered it 
is assumed that cos a - 1, sin a - a and products of 

y y y 

perturbations are negligible in comparison to perturbation 
quantities alone. Under these assumptions equa-tions (8) 
reduce to 


v ip 

and 

Ulp 


1 t 

[Wj + cos a + + U 21 ] sin a 


t t 

[Wi + U 11 ] sin a - [W^cc^ + u 21^ cos 0 ~ a t 


( 9a ) 


(9b ) 


For an oblique shock the tangential velocity component 
is unaltered by the shock whereas the normal velocity com- 
ponent experiences a normal shock. The normal shock rela- 
tion is (reference 8) 


Y+l - z 

Ui P 2 U 1 


U 2 Pi 


v — 1 — ^ 

1 + X-i uj 


= r 


( 10 ) 


or, in differential form 


dU 2 

UJ 


dU 1 


+ 2 


Y-l 1 
r ~~rr- - 1 
Y + l 


dU 2 


U 


( 11 ) 


Now, by definition, U]^ = U^/ai and therefore 
dUj. dU 2 da 2 

= ITT “ -ITT 0 
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Substituting this relation into equation (11) gives 

1 da l 


dU 2 

2r Y ~ 1 - 

1 

dUx 

2r Y ^ - 2 

U 2 

Y + l 


Ul 

Y + l 


a l 


Assuming that 

dUi - U lp - Ujl 


and 


dU 2 « U 2p - U 2 = U 2p - £■ U a 


we have 

U 2 p 


U 2 + dU 2 = ~ Ui 


dU. 


1 + 


or 


U 


2,P 


i Ui 
r 1 


1 + 


2r - 1 


2 J 


1 dUi 


Y+l 

V-l 

Y+l 


2r - 2 


«1 

daj 


Again, for an oblique shock 

V 2 p = ^lp 


dU 


Solving for 




we obtain 


or 


dU x 

~ 

dU i 


U 


2 P 


- 1 


T t 

U 1 1 u 2 1 

- cot C 


w 


Wi 


at 

cot a a - esc a — — 

y W 1 


da j 

Solving for — from first order acoustics we. -have 
a l 

t » 

Pi M 

PI aj 


( 12 ) 


(13a) 


(13b) 


( 14 ) 



and 


For 


1 t ! 


Ti 

Y-l Pl 

= v — = (y - 1) 


Tl 

Y PI 

a i 

the 

assumed perfect 

gas 

2 



a l 

= yRT } 



and therefore 


da ■ 


dT- 


2 T- 


If it is assumed that 
da i 


1^1 
a l 2 


I 

“I 


a l 


dT' 


Ti 


Tl 


, then 
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Now the acoustic particle velocity in Region 1 (upstream of 

I 

the shock front), uj., is related to the perturbation velocity 
components in Region 1 in the following way (see Figure 5) 


i i 

U 11 = U 1 s in (a+0) 

and 

» i 

u 2 i = -Uj cos (a+B). 


(15a) 


(15b) 


Substituting equations (15) Into equation (14) yields 
dUi 


U1 at 

— n — = t; — esc cr cos 0 - a cot a — esc a. 

U i V? i y W i 

The perturbation velocity components in Region 2 at 
the shock front are given by (see Figure 4) 


(16 ) 


= U 2 p cos (<j>+cty) + V 2 p sin (^+a y )«- W 2 


"12 
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and 


i 

u 2 2 = V2p cos ( 4 +ct ) - U 2p sin (4 + a ) 


Noting that is small (by assumption) and that 


W 2 = Wi cos a esc 4 
we have 

t 


ui2= ^2p [cos 4 - sin 43 + 

V2p [sin 4 + cos 43 - 

W i' cos a esc 4 


and 


i _ 

u 22 = V 2 p [cos 4 ~ a y. sin 4>3 - 
U 2p [sin 4 + a y cos ‘f ’ 3 . 

Written in terms of the perturbation quantities in Region 
and neglecting products of perturbations , these become 


and 


u i2 = A 2 u x + B 2 W 1 ot y + C2C1, 


u 2 2 = a 3^i + B 3 W x a + C 3 a , 


where, by definition. 


1 “ 


2r - 1 

Y+l 


o p . 1 »• ■ — 1 

y+l 


2r 1 ~ v - 1 
y+l 


esc a cos 0 + 


cot a 


esc 0 


Y ~ 


2 <n 

, ( y - 1 ) „ 

! ~ r “ 7 + 1 ” 


1 

( 17 a) 
( 17 b ) 
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Ai 

A 2 = sin 0 sin <j> + ™ sin a cos $ 

Bl 

B 2 - sin a sin 4 * + sin a cos 4 
Ci 

C 2 = — sin a cos <j> 
c r 


A3 = sin 0 cos <f> - Aj cos a cos 4 


B 3 = sin a cos 4 - cos a esc 4 > - Bj cos a cos 4 


C 3 = - Cq cos a cos 4 > 


To relate the perturbed pressure in Region 2 at the 

1 

shock to the perturbation properties u 3 , and , the 

pressure rise across an oblique shock front, as given by 
reference 8 , namely 

= _ 2 y - 2 _ y -1 
Pl ytl 1 r+1 J 


( 18 ) 


will be used. In differential form, equation ( 18 ) becomes 

_ 2 

dp 2 dp! 4 Y 0 1 dU x 

?2 Pi 2 yUi-(y- 1) Ul 


As before it is assumed that 
therefore , 

1 

P 2 A4 , C4 

— = — u } + B 4 o + rr- a. 

1 y Wit 

P2 83 


dpi Pi 


PI 


Pi 


and 


dp 2 P 2 


P 2 


P 2 


and , 


( 19 ) 


where 


,-,2 


2yUj - Y ( Y~1 ) 


4y Of 


2yUj - (y-D 


-2 

2yU 


cos 9 , 


A 4 = 


(y-1) 
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B 4 = - 


2 

4y U d 

„ 2 - 

2yU i - (y- 1 ) 


cot <7 


C 4 = 


4y Ui 

2yU j - ( Y -l) 


esc a 


To relate the perturbed density in Region 2 at the 

t 

shock to the perturbation properties uj , a and a the 

y t 

density ratio across the shock, as given by reference 8, 
namely 


P 2 U 1 

Pi " ui 


(20) 


— ' 

P2+P2 Ui+dUi 

- r = U 2 + dUo 

P 1 + P 1 

will be used. Neglecting products of perturbation quantities, 

it can be shown that 

1 1 

P 2 Pi du l dU 2 

P 2 Pi Ul ° 2 


or, substituting for P^/Pl» dd l/ d i and dU2/U2 previously 


derived expressions in terms of uj , a and , 
! 

P 2 a 5 t c 5 


P 2 a l 


u 2 + B 5 ct + 77- o 
y Wp t 


( 21 ) 


where 


o esc a cos 9 Y-l y-l , 

Mj 2 y +1 " 5 


BS = ~ 2 rI^-1 



and 
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c 5 




esc a , 


With the boundary conditions given by equations (17 ), 
(19) and (21) it is now possible to solve for the perturba- 
tion properties in Region 2 in terms of those for Region 1. 


Sound Field in Region 1 . For the plane sound field 
shown in Figure 1 the acoustic particle velocity and acoustic' 
pressure are, respectively, 

WlZl 


U 1 = cos 


a l 


to 1 t 


and 


? 

Pi 


U, 


a l 


toizi 


Pi cos 


- tDl t 


( 22 ) 


(23) 


L a l 


From Figure 5 it is seen that the z^-axis is in recti- 
linear motion in the ^-direction and it lies at an angle of 
(a + 0 - ir/2) with respect to the 5^-axis. Thus, since the 
z j-axis mov'es with speed , 

Z 1 = CCi-Wit] sin (a+6) - m cos (a+0) (24) 


and, therefore, equations (22) and (23) may be written as 

. [ “i 

- si 


Ul = Uq COS 


m cos 


and 


u 0. 

t 1 _ 

pi = y r — pi cos 

a i 


— T [Cl-Wjtd sin (a+0) - 
(a + 0)J - to}t 

r 

— CSl-Witd sin (a+0) - 
L a i L 

m cos (a+0)j - wit 


(25) 


( 26 ) 
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If the transformation x 2 = Si , x 2 = m and U 2 = V7j is 
made, th.en equations (25) and (26) exactly satisfy equation 
(6), Thus equations (3a), (3b) and (6) are consistent with 

linear acoustic theory, 

Sound Field in Region 2 .. For a restricted range -of 
sound wave incidence angles the pressure perturbation down- 
stream of the shock front will be entirely due to acoustic 
waves. This range of angles will now be determined. 

Consider a fluid particle at the point Pj in Figure 6 
at time t. At this same point and time a sound wave inter- 
acts with the shock front and a cylindrical sound wave is 
assumed to be generated. At time t + 5t the plane sound 
wave front in Region 1 intersects (interacts with) the shock 
front at point P 2 and the cylindrical wave has spread a 

radial distance a 2 St from an origin tied to the fluid par- 

j 

tide that has moved from Pj to P 2 , a distance W 2 6t. 

If P 2 lies outside the cylindrical wave front, a re- 
fracted sound wave is generated; if not, a pressure 'wave is 
formed and the original assumption is invalid. Two limiting 
cases exist for 6 for the existence of acoustic waves in 
Region 2 ; namely the two values of 9, say 0^, for which P 2 
lies at the intersection of the cylindrical wave front and 
the shock front. This is shown in Figure 7 for one of the 
critical angles. 

From Figure 7 it can be shown that 

L = [1 + Mi sin (c+0 )] esc 9 a^ 6t 

c c 
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and 

L 


■* 

r 2 

2 

1 i / 2 1 


M2 sin <}> ± 

1 - M2 COS 

* 

] ; 

a 2 fit 


where, in the second equation for L, the (+) is used if 

O<0 <ir and the (-) is used if -ir<8 <0, Equating these 
c c 

expressions results in 


[ 1 + Mi sin (at 0 )3 


esc 0 

c 



[i 


2 2 
M 2 cos 


sin 4 ± 



a 2 

a l 


(27) 


Assuming that 0 is such that acoustic waves do exist 
in Region 2, expressions for the acoustic velocity and pres- 
sure will now be determined. A typical wave front pattern 
is shown in Figure 8. Letting the z 2 -axis be in rectilinear 
motion in the § 2 - cli- rec ' ,: i on with the speed W 2 > the acoustic 
waves are given by 


u 2 


and 

t 

P 2 


U 0 cos 
2 


CO' 


L a 2 


z 2 ~ <o 2 "t 


% . 

Y P2 cos 

32 


CU 2 


a 2 


z 2 - to 2 1 


(28) 


(29) 


Now z 2 may be written in terms of £ 29 02 > t and V7 2 . as 


T t 

z 2 = C?2 ~ W 2 tJ cos (<p-<p ) - r \2 sin (d~d ). (30) 


At the shock front the arguments in Regions 1 and 2 must 
match since the functional form of the acoustic waves was 
assumed to be the same in both of these regions. Thus 
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In summary, the angle of the refracted sound ftave, the 
shift in frequency, and the angles of incidence that yield 
refracted sound waves have been determined. The only remain- 
ing quantity of interest is U 0 > the amplitude of the acoustic 
particle velocity. This can be found only after the shock 
front displacement is determined. 
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Shock Front Displacement, <x(y,t) . The shock front dis- 
placement should be a periodic function of distance along the 
shock front, y, and time, t. From equation (31) the wave- 
length is 


X 


sf 


sin 0 


^ 2 

t 

sin <f> 


(33) 


When acoustic waves exist in Region 2 it is assumed that 


a (y ,t ) = Eq sin 


“1 
L a 2 


= -EjCl + Mq sin (at6)J 


and 


a = E i 

y 




+ M 

0)1 


' ui 


sin 0 cos 

— — 

»1 


- a l 



+ M 


loci ty , 


and 

the sound 


same functi 

onal form. 

Thus 

y - 

[ 1 + 

Ml 

sin 

( at 0 ) 1 
J 

u) L t 

3 






- 

1(34) 



r “1 





“1 

cos 



sin 

6 y - 





L a i 





sin 

(at 8 ) J o> 

1^ 



(35) 

; in 0 

y - 






sin 

( a + G ) j to 

it 

• 


(36) 

i Region 

_2 . 

Since the 

entropy 


in Region 1 is constant (by assumption) and the strength of 
the shock front is variable, the entropy in Region 2 will be 
variable. Thus entropy-vorticity waves are convected down- 
stream with the velocity W 2 . To an observer moving with 


the flow in Region 2, the entropy is a function of position 
alone since the acoustic compression-rarefaction process has 
been assumed to be isentropic. 
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The entropy- vorti city wave front in Region 2 generated 

1 

J 

by a plane sound wave in Region 1 is shown in Figure 9. The 
E-V wave front ray. angle, 3, is given by 

Mj. cos a cot <f> sin 9 


B = t an 


-1 


j 1 + Mi sin o cos 0 J 


(37) 


The E-V wavelength is 


X = Xl 5424 

evw - 1 sm 9 


(38) 


and since the Z 3 -axis moves rectilinearly in the 5 2 -c lir ec tion 


z 3 = ” ^2*3 cos (<p-3) - T)2 sin ($-3) . 


(39) 


Thus the argument of the entropy-vorticity wave, 

2irzq/A , can be shown to be 
° evw’ 

2 ttz 3 27rsin8^[C2~ w 2' t ] cos ^ , } , "0) -r l2 s i n ( < S , ~&)J 


evw 


A sin 3 


(40) 


Perturbation Velocity Components in Region 2 . The per- 
turbation velocity components in Region 2 consist of two 
pants: one due to the sound field and one due to the vofti- 

city waves. Assuming the vorticity waves to be of' the same 
functional form as the sound waves, the perturbation velocity 
components may be written as 


U32 


U o cos 
2 


0)2 


L ^2 


w 2 t 


i 

cos (<£>-<(> ) + 


F j cos 


2ir z 3 " 
~ 

evw J 


(41) 



and 


i 

u 22 



as 


(42) 


Perturbation Pressure in Region 2 , From equations (3a) 
and (3b) it can be shown that the perturbation pressure in 
Region 2 is only that associated with the sound field; that 
is 

' Uo 

P 2 2 w 2 

= y cos , — z 2 - w 2 t . (29) 

P 2 32 L 32 


Perturbation Density in Region 2 . With only an acoustic 
pressure perturbation in Region 2 the perturbation density is 
composed of an acoustic part and an ent ropy-vort i c i ty part; 
that is 

P 2 °2 w 2 ’ 21 TZ 3 

■ — = cos z 2 - w 2 t + H; cos — y . (43) 

p 2 3 2 L a 2 J L e vw. 

In review, there are now five unknown quantities in 

Region 2; namely, Ej, Fi, Gj , Hi and Uq . There are only 

2 

four equations relating these to known quantities at this 
point in the analysis and, therefore, an additional relation, 
the continuity equation in Region 2, will now be derived. 

Continuity Equation in Region 2 . Consider a change in 
variables in Region 2 such that 


?3 ~ £2 ~ W 2 t 
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and 

n 3 = n 2 j 


that is , a change to a coordinate system that moves with the 
unperturbed flow. Now, in terms of £3 and q 3 , 

t 1 

z 2 = £3 cos - ri 3 sin ) 

and 

z 3 = £3 cos ( 4 >- B ) “ n 3 sin ($-B). 


Equation (2), written in dimensionless form for the £ 3 , n 3 
coordinate system, is 

ap 2 u 2 1 3 p 2 2 1 ° p 2 

~ „ — cos ($_(}> ) — sin C ) — c — + 

3t a 2 3 £ 3 a 2 3n 3 

16 116 11 


U 0 


P 2 + 


2 ~ i 

- P 2 


a 2 


^ I ~ t 

» 9u 12 i 3u 22 

cos ) — — — - sin (<|>-<j> ) 


S 0 . 


u 3 


3 n 3 J 


16 1 11 11 


Retaining only those terms that are of the order of magnitude 
of one, this becomes, in dimensional form, 

3 p 2 


St 


+ P 2 


1 t 

SUj2. 3 u 2 2 

3 £ 3 3 n 3 


= 0 , 


(44) 


T t < 

Substituting the expressions for p 2 , u 32 and u 22 into equa- 
tion (44) yields 


G x = Fi cot (<j>~0 ) . 


(45) 



There are now four unknown quantities in Region 2 and 
four equations thereby forming a determinable system. The 
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unknowns will now be determined. 


Evaluation of Ey , Fj , Hi and Uq . The following equa- 

2 

tions are to be solved simultaneously at the shock front: 


t ! 

A 2 uj + B 2 W i + C 2 a t = u i2 , 


i i 

A 3 u x + B 3 W x a y + C 3 a t = u 22 , 


i P 2 

Ul t Bk a + rr— = 

1 V Wit 

*1 P2 


A 5 , C 5 p 2 

— *- ui + Be a + 77— a, = 

1 0 y Wi t 

a i 1 P2 


It is also known that 


■ “>2 1 T 

Uq cos — z 2 - m 2 t cos ) + 

2 L a 2 


F i cos 


2tt z 3 


u 22 = - Uq COS 

2 


z 2 ~ w 2 t sin C4~4> ) + 


' 2m z 3 
Gi cos y — ■» — - 


P 2 2 r w 2 ' 

=■ y cos — — z 2 - W 2 t , 

P 2 a 2 L a 2 



T 

P 2 

u ° 2 r 

W2 


- 



" 2ir z 3 _ 


— „ ^ „ 


z 2 - 

- 2 t 


cos 


_ 

“• L- y 0 

_ 


\ 

p 2 

a 2 

- 32 


" 

n 


L eyw J 

t 


-&)1 






Ul 

= Un COS 

O 

— — 

21 

- Ci>l t 

3 




4 

-Hi 


- 





“1 



r “1 




a 

= Ei • — sin 0 

COS 

sin 0 y 

- 


y 

31 



l 31 







Mi 

sin (cr+0) 

Wit 

9 


a = - Ei [1 + Mi sin (ai-6)] toi cos 


“1 


L ai 


sin 0 y - 


£l + Mi sin (a+6)j oiit 


and 


Gi = Fi cot ($-$) . 


31 

(43) 


(22 ) 


(36) 


(35) 


(45) 


It can be shown that the arguments of the cosine terms 
in equations (22), (29), (35), (36), (41), (42) and (43) are 
all equal at the shock front. Thus equations (17a), (17'b), 

(19) and (21) may be written as 

“ 1 

A 2 U 0 + B 2 Wi Ei — sin O — C 2 wi Ei Cl + 


Mi sin (crt0)] = Uq cos (^>-<j> ) + F^, (46) 

2 

w 1 

A 3 U Q + B 3 W . E 1 — sin 0 - C 3 m Ei Cl + 

1 a , 


M, sin (ate)] = - U 0 sin ( <}>— ) + Fi cot ($-0) 

2 


( 47 ) 


5 



32 


0 2 <u 1 c Lj 

A4 + Bj| Ei — — sin B - 77— wj Ei [1 + 

Si Sj 1 


Mi sin (0+9)3 - y 


1 ’“1 c 5 

A5 - — t. — + B5 Ei ■ ■■ ■ sin 9 — 77 — <1)2 Ei [1 + 
ai ai 1 


Mi sin (o+ 0)3 = • — — + Hi 

32 

These four equations can be arranged to give 


An E 1 + A I2 f 1 + A 13 U 0 + a 14 Hi = Dj 

2 


A 2 1 El + A 2 2 F 3 + A 2 3 u 0 + A 2 4 Hi = D 2 

2 


A 3 1 Ei + A32 Fi + A33 U 0 + A34 Hi = D3 

2 


A41 Ei + Aif .2 F 1 + A43 Uq + A44 Hi = D4 

2 


11 * 


1 + Mi sin (a+6) j C 2 - Mi sin 0 B 2 wi 


A l 2 " 1 


A 1 3 = cos ($-<!> ) 


A 14 = 0 


1 + Mi sin (0+0) J C3 - Mi sin 0 B3 


A 2 2 = cot ( 4> — 3 ) 
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A 23 

a 24 

a 3 i 

a 32 
A 3 3 
A-34 
A 4 1 
A 4 2 
A 4 3 
A tf 4 
Di 
D 2 


T 

- sin (4>-<j> ) 


0 


+ Mi sin 


(o+e) ] 


C 4 - Mi sin 9 B 4 j 


«1 


0 

aj 

Y ~ Mi 
a 2 

0 


+ Mi sin 


( cr+ 9 ) j C5 


Mi sin 9 B5 


“1 


0 

a 1 

a 2 
1 


A 2 U0 1 

A 3 OOj 


D3 - A 4 Mi U 0 
D 4 “ A 5 M 1 u 0 


Computations . To compute the perturbed flow field in 
Region 2 due to the interaction of plane sound waves from 
Region 1 with an oblique shock front, values for the follow- 
ing quantities must be known: 

_ _ t 

ai ,y,Mi ,0 , 6 ,r ,a 2 ,<f> , 9 ,U 0 ^ ,u)i ,<J> and 
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To illustrate the computational method, y, aj (or Tj 
a, 0, Uq and Wj will be treated as independent variables, 
by-step procedure is as follows 


1. Select a temperature Tj. 

2. Select a working fluid, i.e. y; e.g., air with y = 1.4 

3. Select a Mach number Mj. 

4. Select an oblique shock angle a. 

5. Select a sound wave incidence angle 0 (if possible, i. 

consider 0 ’s). 
c 

6. Select a sound wave particle velocity U Q . 

1 

7. Select a sound wave circular frequency toj. 

8 . Compute r as follows (reference 8 ) 


’l P 


2 (y+1) M 2 sin 2 a 


2 P 1 2^1+ M 2 sin 2 a J 


9. Compute d as follows (reference 8). 


& - a - tan" 


1 . I 

— tan a 
. r 


10. Compute as follows (reference 8): 


VyRt" 


V t 2 /t 1 


and 


t 2 /Ti 


[ 1 + 1 r 


(y+1) 2 

2(y-l) 


TP- 


, Mj, 

The step- 


. must first 


(51) 


(52) 


(53) 
(54 ) 
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11 . 


12 . 


13. 


1M-. 


15 0 


where 

Uj = Mj sin a 


C 5 5 ) 


Compute <J> by noting from Figure 1 that 


<j> = ^ + 6 - a 


C 5 6 ) 


Compute 0^ from equation (27), i.e. 


El + Mi sin (a+0 )U esc 0 

c c 


Mi cos a + 


r “ 2 M 2 . 2 “1 1/2 

a2 Mi sm a ' 


-2 

a l 


where (+) is used if 0 < 0 < ir is being considered 
or (-) if 0 > 0 > - ir is being considered. (This is 
a transcendental equation for 0^ and an iterative method 
of solution is required.) 

Check to see if the chosen value for 0 is permissible, 
i.e. | 0 | < 1 0 c | . 

Compute 0 from equation (37), i.e. 


g = tan 


-1 


Mi cos a cot <j> sin 0 
1 + Mi sin a cos 0 


Compute <p from equation (32), i.e, 


a 2 sin 0 
ai sin $ 


1 t Mi sin (at 9) 

. - 

1 t M 2 cos ( <J> — tj> ) 


where 


cos a 
sin <J> 


a i 


a 2 


M 2 = Mi 


(S'7) 
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(This is a transcendental equation for and an itera- 
tive method of solution is required. ) 

16 0 Compute A.. and D. 0 (At this point the A..*s and D.’s 
lj 1 1 

can be computed using their defining equations.) 

17. Compute Ei as follows: 

a. Note that each D. contains Un and that Mi will be 

1 1 

a factor of the coefficient determinant (solution 
by Cramer's rule). Therefore it is convenient to 
define new D^'s and concomitantly a new Ei; i.e. 


D . = D . U 0 
i i 


and 


Ei = Ei Ug^/mi . 

b. Using Cramer's rule the solution of equations (50) 
for Ei is 


Ei = [D 1 A 22 A 33 + E 2 A 32 A 13 + D 3 A 2 3 A • ’ - 

12 

A A A 

D1A23A32 - E 2A 1 2 A33 - D3A22A133/EET [A^J 


(58) 


where 


DET [A„] ~ A 11 A 22 A 33 + A 21 A 32 A 13 + A 31 A 23 A 12 ~ 

A 31 A 22 A 13 “ A 2 1 A i^A 3 3 - AiiA23 A 32 (59 ) 


and finally 


El 


Ei Uq^/wi . 


18 0 Compute Fi as follows: 
a. Define a new fi as 
E 1 = ?* U 0i . 
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19 . 


20 . 


b. Using Cramer's rule the solution of equations (50) 
for is 

Fl = [D1A23A31 + D 2 A 11 A 3a + ^ 3 a 13 a 21 “ 


D1A21A33 - D 2 A 3 1 A j 3 - D 3 A 1 1 A 23 ] /DET [A..], 

(60) 

Therefore 

Fi = r* u 0 . . 


(It should be noted that is a factor of both the 
numerator and denominator and, therefore, Fj is 
independent of frequency.) 

Compute U 0 as follows: 

a. Define a new Uq as 

2 

U 0 = Uo U 0 , • 

2 2 1 

b. Using Cramer's rule the solution of equations (50) 

A 

for Uo is 
2 

sV it ti 5V sV 

Uo = [DiA 32 A 2 i + D 2 Ai 2 A 3 i + D 3 AiiA 22 - FiA 3 iA 22 - 
2 


D 2 A l l A 3 2 - D 3 A 2 iAi 2 !]/DET C A ^ H 


(61) 


and finally 


U 0 = U 0 U 0 o 
2 2 1 

(Note that Uq is also independent of frequency.) 
2 

Compute Hi as follows: 
a. Define a new Hi as 
Hi = Hi Uq o 
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21. 


b 0 Using Cramer's rule the solution of equations (50) 
for Hi is 

A C J. 

H 1 = l[ A 2^A32A43 + A 3 i A 4 2; A 2 3 + A41A33A22 ~ 

A 4 1 A 3 2A 2 3 - A 31 A 22 Ai f 3 . - A21A33A42] + 

D2 [a 1 1A32A4 3 + A12A33A41 + A 13 A42A 3 1 - 

A41A32A13 - A31A12A43 - AiiA 33 A4 2 j - 

.1, 

D 3 [AhA 2 2 A 43 + A12A23A41 + A 3 iA 4 2 A 2 i - 
A41A22A13 ~ A 2 iA 12 A 4 3 - AHA23A42] + 


D 4 DET [A . .1 

and, therefore, 

Hi - Hi U 0 . 


/DET [A . . ] 
3-1 


( 62 ) 


(Note that Hi is independent of frequency.) 

T ! 

Compute the acoustic pressure ratio p 2 /pi as follows: 
a. From equations (23) and (29) 

w 2 


1 u o P 2 a l cos 
P2 2 


z 2 - u 2 t 


a 2 


T 

PI 


Dfl Pi a 2 GOS 


®1 


a l 


Zi - Wit 


or upon taking the root-t ime-mean-square of p 2 


end pi 


Cp 2 l 


U, 


rms 


2 P 2 a 1 


Cpi3 


u 0 


( 63 ) 


rms 


1 P 1 a 2 
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where 


U 0 


2 




5 


P 2 
PI 


= 2 V 

y+ 1 


2 

Mi 


. ^ 

sm a 


y-1 / ~ 

- — =- reference 
yt 1 


8), 


and 

a 1 

— — is as given in Step 10 . 
a 2 


b. From the definition of sound pressure level, the 
change in SPL across the shock is 


SPI >2 - SPLi = 20 log i q 


CP23 


rms 


Cpi3 


rms 


and therefore 


SPL 2 - SPLi = 20 logig 




2 P 2 a 1 


U« 


1 Pi a 2 


22. Compute the ratio of the acoustic particle velocitie 

t I 

U 2 /U 1 from equations (22) and (28); i,e. 

o > 2 


, Uq cos 
u 2 2 


a 2 - m 2 t 


a 2 


1 

“1 


U 0 cos 
1 


0)1 

- — Z 1 - 0)1 t 

31 


I 1 

or upon taking the root-mean-square of Ui and u 2 


[u 2 ] 

[uj] 


rms 


rms 


U 0 



2 



% 

= 

u„ 2 


(64) 


(65) 


( 66 ) 
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23. Compute the intensity of -the vorticity (shear) wave , 


[u ] /W?, as follows: 

evw z ’ 

rms 


a. From equations (41) and (42), the velocity pertur- 
bation due to the E-V wave is 


' r 2 

u = [Fi t GiJ 
evw 1 1 


1/2 I" 2ir z 3 


cos 


evw 


or taking the root-mean-square of u 


evw 


[u ] 
evw 


rms 


2 2 1 X / 2 
F i + G i 


2 2 2 

Now G^ = F i cot (<j>-3) and therefore 


2 2 2 r 2 , 

Fj + Gj = ?i [1 + cot (<j>-3)] 


[u ] = Fj esc (a -$)// 2 


evw 


rms 


U, 


= Fi esc (o-3) 


/2 


c. From equation (57) 


Wo = Mi a 


cos a 


1 ai — : X 

1 1 sm <p 


and therefore 


Eu ] 

evw 


M/ 


rms 


W* 


* esc ( g- g ) sin <j> 

/2 Mi COS * 1 


(67 


24. Compute the perturbation entropy in Region. 2 as follows 
a. For a perfect gas 


ds = c ^2- - c — 


v p p -p 



25 . 


26 . 
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Replacing dif f erential quantities by perturbations, 
this becomes 


t 


t 

P 2 


i 

P 2 


S? “ C — - c — 
^ V - p - 

P 2 P 2 


( 68 ) 


b. Substituting equations (29) and (43) for P 2 /P 2 an ^ 

t _ 

P 2 /P 2 s respectively, yields 

U<> 


S2 - y c 


cos 


v - 
a 2 

Uo 


c cos 

p a 2 


c Hi cos 
P 


w 2 


z 2 - 0)2^ 


a 2 

W 2 

L a 2 

2 ir z 3 
~X 

evw 


z 2 “ 


or after noting that y c = c 

v p 


t 

s 2 * 

= - Hi U n cos 

c 1 
P 


2 IT z 3 


evw 


(69) 


Compute the perturbation density in Region 2 from 
equation (43); i.e. 


1 




cos 


a 2 


U) 2 

z 2 ~ u 2 1 

a 2 


+ H 3 cos 


' 2 tt z 3 ~ 


X 


L evw _ 



U, 


(70) 


Compute the shock front perturbation velocity, , 
from equation (35); i.e. 

sin 8 y - 

. a l 

£l + Hx sin (cr-f- 6 )J w * 1 1 (71) 


“t 

U 0 . 


= - EiCl + sin(a+ 0 )j cos 
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or in terms of the ratio of the root -mean-squares of a^_ and 


Ui 


01 

t rms 


[Uj] 


Ej ) jl + H 1 sin (cf+0 ) 


rms 


27. Compute the peak-to-peak amplitude of the shock front displace- 
ment , 1 2a J , from equation ( 34 ) , i . e . 


| 2a | = 2E 3 = 2E" U 0 /u, 

1 


or 

I 2cc I 03 1 


A FORTRAN IV program to make these computations is given in 
Appendix A. 

Results 

Rather than making computations for arbitrary shock fronts, it 
was decided to illustrate the sound wave - shock front interaction 
process for two Saturn V interstage flare angles, namely, 5-8° 58' 
corresponding to the flare between the Service Module and the S-IVB 
stage, and 6 = 16° 44', corresponding to the flare between the S-IVB 
and the S-II stage. 



( 72 ) ' 


(73)’ 
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Appendix B contains computed data for the following cases 

M 1 1.50 2.00 2.50 3.00 1.75 2.00 2.50 3.00 

6 8.98 8.96 8.97 8.98 16.74 16.73 16.74 16.75 

The data column headings are defined as • 

THETA - 0 , the sound wave incidence angle , in degrees 
BETA - 8, the entropy-vorticity wave angle, in degrees 
PHIP - <j> , the refracted sound wave angle , in degrees 

1 I 

PP2/PP1 - P 2 /p 1 , the ratio of the acoustic pressure in Region 2 to 
that in Region 1 evaluated at the shock front 

SPL2/1 - SPL^ - SPLj, the change in sound pressure level across the 
shock, m dB 

U02/U01 - U Q /U 0 
2 1 

UEVW2 - [u' 3 /(M Q W ), vorticity wave intensity 

evw rms 1 

.SPRCU - [s'] /c U n 
2 shock P 

J 

AT/UP1 - [a /uj] , the ratio of the shock front velocity to the sound 

shock 

particle velocity in Region 1, evaluated at the shock front 

f 

PPAMPL - [2ooj 1 /u 1 ] , the ratio of the product of twice the shock front 

shock 

displacement and the circular frequency in Region 1 to the sound 
particle velocity in Region 1, evaluated at the shock front. 

It should be noted that shock front resonance is predicted for all 


but one of the cases considered, namely, M 1 = 1.50 and 6 = 8.98: The data 

are not displayed graphically . 
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Bp Shear Wave ~ Shock Front Interaction 

Since shear waves (vorticity waves ) are generated when 
plane sound waves interact with an oblique shock front and 
since a given flow often interacts with more than one oblique 
shock, a method to determine the flow field produced by shear 
wave — shock front interaction is required. Rather than de- 
velop a solution for an oblique shock, a coordinate transfor- 
mation will be performed so that Ribner's [1] analysis of ' 
the shear wave — normal shock front interaction can be used. 

Mathematical Formulation of the Problem 

Unsteady-Flow Problem ., Consider a uniform velocity field 

W A perturbed by the addition of a sinusoidal shear wave of 

arbitrary amplitude, wavelength and orientation. This flow 

undergoes an oblique shock as shown in Figure 10 . Since , in 

general the perturbation velocity w is not parallel to W. 

and since the s'hear wave pattern is frozen in the uniform 

flow, the perturbed flow will be unsteady from the viewpoint 

of a stationary observer on the shock front. For example, 

the nodal lines N-N will be convected into the front and 

s 

the nodal points N g will move up the front for the case de- 
picted in Figure 10. 

Equivalent Steady-Flow Problem . If an observer moves up 

the shock front with the speed of the nodal points, , 

s 

flow will appear to be steady. To determine this speed 


the 



45 


consider the nodal line Mj-H at time t as shown in Figure 
11. This line is convected into the shock front by the uni- 
form flow W . At time t + fit the point has reached the 

A 

* 

shock and now represents the nodal point N or N to dis- 

S 1 . S! 

tinguish it from the nodal point at time t. 

From Figure 11 it can be shown that 


NtN = NhN cos [90-(6 + c) 3 + N N cos o, 
A s i s l s l s l 


NiN sin [9O-(0+a)3 = N N sin a, 
si si si 


NiN 1 S W A fit and N N* = V XT fit 
si A sisi 


Solving these equations for yields 


V„ = W. cos a - W. sm a tan 8 
N A A 

s 


(74) 


Moving up the shock front with the speed the flow 

s 

appears to be steady and the unperturbed flow appears to have 

the velocity W as shown in Figure 12. The angle that W 
A A 


makes with the normal to th'e shock front is 
-1 


9 A = tan 


sin a tan 


. W. sin a 
L A 


or 


e . = o . 

A 


Thus from this viewpoint the unperturbed uniform flow 


is aligned with the perturbation velocity component « . The 

"A 
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magnitude of is given by 

W = W sin o sec 0 (75) 

A A 

and therefore it gets very large as 8 approaches 90 degrees. 

The flow field is now in Ribner's equivalent steady- 
flow form. His solution will be stated next and it will then 
be transformed back to the original oblique shock form. 

Ribner's Solution 

Ribner’s coordinate systems and velocity components are 
shown in Figure 13, In terms of the oblique shock proper- 
ties, his velocity components before the shock front are 

W . — W , sin o sec 0, (76) 

A A 

U A = W sin a (77) 

and 

V = W sin a tan 0, (78) 

A 

His perturbation velocity component w^ is 
w A = e W A cos (kn A ) 

or in terms of oblique shock quantities 

w A = s W A sin o sec 0 (79) 

From this last relation it should be noted that the 
assumed amplitude of the perturbation is strongly dependent 
on 0; i.e. it goes as sec 0 „ Thus for a given W^ } 


a and e 
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(to be taken as a number less than about 0.1) the strength 
of the perturbation varies with 0 and it is unbounded at 
0 = 90 degrees. This point will be discussed further when 
the flow field is transformed back to its original unsteady 
form . 


The Critical Incidence Angle . Two forms of the solution 
for all flow quantities appear, one for W < 1 and the other 
for W > 1. W = 1 forms the dividing line and, since W depends 
on both U. and the inclination 0, the sonic conditions gives 
a critical value of 6 in terms of 0^. From Ribner [l] 


6 = ± 

cr Xl 


2m 


(80 ) 


where 


U, 


m = 


ni o 

2 A 
v-1 —2 

1 + V °A 


, U A “A . 

U. = = — sm a , 

A a A a A 


Thus for a given W and a there exists a 0 


cr 


Velocity Perturbation Field, W < 1. From Ribner [1] 


w 


w 


= S cos [k cos 0 (y-x tan <f> ) + 6 3 + 
** s 


A 


n cos [k cos 9 (y-x tan <j> ) + 6^] 


(81a) 


The bar superscript 
ratio of the_veloci 
sound; e.g. W = W/a 


denotes the 
ty component 

» D A = U A /a A 


Mach number formed by the 
and the local speed of 
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w 


W , 


-, = 8 n sin [k cos 8 (y-x tan <£) + <$] 
w p 


(81b) 


where 


w . | 5 W. t cos 0 = e W, sin o 
A 1 A A 


* 5 5 V 


2 2 
A + B 


/~ 2 ~ 2 2 

R H ~ —2 — -1 — — exp [- x k cos 8 8/8] 
m p w 


6 = tan 

s 


-1 


-B 

A 


-1 


<f> = tan [m tan 0] 


tan 


-1 


r _ 2 

U t an <j> 
2 
e 


6 = t an 

P 


-1 


c 8 -d tan 
w 


L d 8 +c tan i> 
w 


2 _2 
8 = 1 - U 


2 _ 2 

8 = 1 - W 

w 


~ ( m— 1 ) 

A = sec $ + 2 (m-1) cos d> + a • — — sin <j> 


m 


(m-1 ) 


~ D 

m 


= m 

CE 

+ 

DF 

2 


2 


C 

+ 

D 

= m 

CF 

_ 

DE 

2 


2 


C 

+ 

D 


sin <j> 




Ribner gives the denominator as dB - c tan 4> but a positive 
sign was obtained by this author. W 
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c 


d 


C 



i 


- F 



m 


[ 


Y-l 

y+1 


+ 


3-v 

Y+1 


m 


j tan $ 
2 (m-1 ) I 

yn J 


- £ (m-1 ) + 

sin $ cos <}> 


D 


E 


F 


8 2 Pj. r 

— ~ (m-1) [1 + (m-1) cos 4>] = — -r D 

2 6 


3 


2 [ 1 - ?TT" 


| + 2 (m-1) — 


2 2 . 

8 cos <j> 


[2 ( m-1 ) sxn (j> cos <|> ] t 


w 


Velocity Perturbation Field, w >- 1. From Ribner [1] 

A A 

■ A = S cos [k cos 6 (y-x tan 4> ) ] + 

% A t 

n"' cos [k cos 0 (y-x tan $ )3 

t f 1 

= 8 ]f'"bos[k cos 6 (y-x tan 4 )3 
w 

where 

|w„ I" = W„ e cos 6 = e W sin a 
1 A 1 A A 





S = A/m 


(82a) 

(82b) 
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_ c_ s m y 

^ IB f x. \ " 

cos C (f> — p ) 


(Ji = | [ | tj) j - p 3 | with sign of 

\i = cot -1 6 I" g e V w i 2 - 1 1 
w L w J 


A = A for W < 1 


a = m 


! f 

C + GF 

t t 

E + GD 


b = 0 


~ = i r 

c = - D - F 
m 


C = 2 


Y-l 

y + 1 


m - 2 [1 + (m-1) cos <j> ] 


i 2 

D = (m-1) [1 + (m-l) cos <f>J 


E = (m-1) sin 4> cos $ - £ 1 + in J tan <j> 


F =2 (m-1) sin $ cos <b 


1-3 tan <|> 


G 5 


w 


3 + tan <j> 

VI 


= tan (y-<j>) 


Pressure Perturbation Field., W < 1. From Ribner [1] it 


can be shown that 

“ e 7vS '' Y ~ T - ITh 11 cos Ck cos 0 (y-x tan 4> ) + 6 ] 

p (Y~l)~ (y+1 ) m J Y p 

___ (83) 

+ Instead of cos (<j>-y) Ribner has sin (<j> + y). This is in error 
since the cos (<j>-y) follows directly from his equations (32) 
and y 5 cot -1 - 3 . 
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where n 3 <J> and 6 are the same as those defined for the 

P 

"subsonic" velocity perturbation case. 

Perturbation Pressure Field, W > 1 . From Ribner [ 1 ] 

6p 2 y is sec if TT " * r . Q , . 

p = 8 TF T?- T Ftn ~ m ' 11 008 [k 008 8 (y -* tan * >] 

(84 ) 

* :V 1 

where n and tj> are the same as those defined for the 
"supersonic" velocity perturbation case. 

Shock - Wave Perturbation, W < 1 . From Ribner [ 1 ], the 
local deflection fix from the plane x = 0 is given by 

6x = 7T~\/ a + b cos [k cos 8 y + 6 , , ] ( 85 ) 

k cos 6 V J shock 

where a, b are the same as those defined for the "subsonic" 

velocity pert ur bat ion case and 

where 

_ 1 ~ ~ 

6 = tan (a/b). 

shock 

Shock - Wave Perturbation, W > 1 . From Ribner [13 

< 5 x = t — - — — — 7- sin [k cos 0 y] (86) 

k cos 8 J 

where a is the same as that defined for the "supersonic" 
velocity perturbation case. 
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Oblique Shock Solution 


The following simple coordinate transformation will take 
the flow from Ribner's equivalent steady flow to that of an 
unsteady flow as viewed by an observer fixed to the shock 
front 


x = x 


and 


or 


y = y + V N t 


X = X 


and 



s 


where, from equation (74), 


(87a) 


(87b) 


V 


N 


s 


[cos a - sin a tan 8]. 


(74) 


Thus x is position measured normal to the unperturbed shock 
front with x = 0 representing the rest position of the shock 
and y is position measured parallel to the shock front rest 
position; y = 0 is an arbitrary position. 

The solution, m terms of the oblique shock properties, 
will now be given. 


The Critical Incidence Angle . For the oblique shock, 
equation (80) is 


0 

cr 


tan 


• 1 / ( yt 1 ) (m- 1 ) 

V 2m 


(80) 
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where 


m = 


v + 1 ~ 2 2 

■^2 — sm ° 

1 + +-J- M a sm a 


W, 


and 


A ~ a, 


6 is measured from the negative x-axis 
cr 


If the critical incidence angle is to be measured from 
the horizontal, then from Figure 14 , 


0 =90-0-6 

cr cr 


( 88 ) 


The Shear - Wave Incidence Angle . If 0 is measured 
from the horizontal instead of the negative x-axis then, as 
shown in Figure 14, 0 becomes 


6 = 90 - a - 6 


( 89 ) 


Velocity Perturbation Field, 0 < 0 < 0 . Upon 

cr. cr c 

& u__ 

transforming to x and y and expressing the velocity pertur- 
bations in terms of streamline coordinates E, and n (see 
Figure 14 ), equations (81) become 


w 


w 


A 


- S sin o cos (4>-6) cos [k cos 0(y-x tan 4>-V^ P ) + 


6 ] + H sin o cos [k cos 0(y-x tan tf> -V^ t) + 


+ 0 =90-o 

Cr * 


0 = tan 

cr 1 


<$ + 4> - <5 ] 

P 


0 and 0 =90-0+0 where 

cr 1 cr 1 cr 1 

u 


(90a) 


-1-, /(y + lMm-l) 


2m 
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= S sm o sin(<f>-6) cos [k cos 0(y-x tan <|>-V^ t) + 

s 

~ ^ f 

6 ] + n sin o sxn [k cos 0(y-x tan © -V„ t) + 
s N 

s 


6 ^ + <}> - 6 ] 


(90b ) 


where 


! W A ' H £ W A 


<f> = 6 + o - 90 


and where the other quantities are as previously defined for 
this case; i.e, W < 1. 


Velocity Perturbation Field, 0 < 6 or 0 > 6 

cr, cr 

£ u_ 

Performing these transformations for the "supersonic" case 
equations (82) become 


S sin a cos (4>-$) cos [k cos 0(y-x tan t) + 

s 

A A ~ ~ t 

5 ] + n sin a cos [k cos 0(y-x tan d> -V,, t) + 
s J T N 

s 


6 +4-6] 

P 


Ola) 


S sin a sm (<j>-<$) cos [k cos 8(y-x tan t) + 

s 

* ~ ~ t 

6 ] + n sin a sin [k cos 0(y-x tan <f> -V t) + 
s N 

S 


6 + $ - 6 ] 
P 


(91b) 
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where 



$ = $ + cr-90 

and where the other quantities are as previously defined for 
this case; i.e, W > 1. 


Pressure Perturbation Field, 


transforming to x and y , equation 



(83) 


< e < e 

cr 

u 


becomes 


Upon 


Sp _ 2 y m sec $ 

p £ ( y - 1 ) - ( y+ 1 ) m 

5 ] 

P 


n 


- ~ I 

cos [k cos 0(y-x tan $ -V^ t) + 

s 

(92 ) 


where H , 
"subsonic" 


$ and 5 are the same as those defined for the 
P 

velocity perturbation case. 


Pressure Perturbation Field, 6 < 0 or 6 > 8 

cr. cr 

- £ u 

After transforming to x and y, equation (84) becomes 


5p 2 y m sec <p ** r . ~ * 

p = e Ty'-I ) - (y'+ 1 ) rn n COS Ck cos e( y~ x tan + - V N 


(93) 


where k and $ are the same as those defined for the 
"supersonic" velocity perturbation case. 


Shock - Wave Perturbation, 0 <■ 0 < 8 

cr. cr 

l u 


Transforming 


x to x, 6 x to 6x and y to y - t, equation (85) becomes 
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Sx - E — — 1 / a 2 + b 2 cos [k cos 0(y-V t) + <5 ] 

k cos Q \ J N shock 

s 


(94) 


where a, b and 5 are those defined for equation (85). 

SilOCK 


Shock - Wave Perturbation, 0 < 6 or 0 > 0 . Up'on trans- 


cr 


cr 


forming to x and y, equation (86) becomes 


Sx - sin [k cos 0(y - V t)] (95) 

k cos 0 N 

s 

where a is the same as that defined for equation (86). 

* 

The most important properties of the perturbed flow field down- 
stream of the oblique shock front have been presented. However, there 
remains one property of interest that is normally viewed by a station- 
ary observer on a surface inclined at the angle 6 - the property is 
pressure and the observer is usually a small transducer which is taken 
herein to be much smaller than the shortest wavelength to be considered. 
Thus , in the next section the pressure at the transducer will be related 
to the shear-wave w^. 


Perturbed Pressure at a Stationary Transducer . The system to be 

treated is shown in Figure 15. It can be shown that the surface wave 

number k is related to the shear - wave wave number k by 
sur 

I 

i- = cos (ottji -6) cos 0 

COS (j) 


sur 


t 


k 


(96) 
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or since the wavelength and wave number are simply related' 

(kA = 2tt = k A ) 
sur sur 

i 


A 


sur 


cos <j> 

t 

cos (a+<j>. -5 ) cos 6 


(97) 


where A is the shear - wave wavelength. 
The period of oscillation is 

P = 2 71 

k c o s 0 V , t 
N 

s 


(98) 


hnd therefore the frequency is 


f = 1/P , v M 


or substituting for V from equation (74) 

s 
W 


f = cos ( o + 9 ) 


(99 ) 


It should be noted that f = 0 when 0 
pressure wave fronts are stationary. 
The surface pressure wave speed 


V 

p 

e sur 


A 


sur 


P 


= 90 - a ; i 


, V , is 
r sur 


e . 


the 


or 


V 

p 

r sur 


t 

cos tp cos (a+0) 

i 

cos 9 cos (a+$ -6) 



(100) 


or, if V is expressed in terms of the local isentropic 

P S ur 

speed of sound, 


V 

p 

r sur 

a 


I 3 

cos <j> cos (a + 9) ~ A 

— 1 T M . — 

cos 9 cos (a+<j> -5) 3 


(101) 
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where 


M 

A 



a 



5 



and 

T 

A 

Y~ is given by equation (54) with T = T , T = T and 

M sin a = U, . 

A 1 

Only the root-mean-square of the pressure remains to be evaluated. 
Taking the rms of equations (92) and (93) yields 


(dp) 


rms 


2 y m sec <j> 

(y+1 )m- (y-1 ) np 


( 102 ) 


where II is either II* or II** depending on the value of 0; i.e. whether 
W < 1 or W > 1. 

The "sound" pressure level is, by definition, 

(dp) 

SPL = 20 log 1Q EEE dB (103) 

P ref 

where p ^ = 2 x 10~ 4 y bar - 0.290 x 10" 8 psi. 
ref 

Computations . To illustrate a computational method, y, a^ (or T^), 

M , a, 0, e and k will be treated as independent variables. Numerical 

A * 

computations using the following steps are given in Appendix D. 

1. Select a temperature T . 

A 

2. Select a working fluid, i.e. y. 

3. Select a Mach number M . 

A 
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4 . Select an oblique shock angle o . 

5 . Select a shear wave incidence angle 0 . 

6 . Select a shear wave amplitude e . 

7. Select a wave number k. 

8. Compute m as follows (reference 8) 


m 



p 



(-Y+1) Ufl 

2 + (v-1) U 2 
A 


where 


(104) 


U = M sin a 
A A 


9. Compute a as follows (reference 8) 


a=a fl /r7T^ 


T_ 

T 


r 2 + (Y-i) u 2 i r 2 y_ if _ 

L A j L v-1 A 


(105) 


(y+i) rr 2 
(y-i) a 


10. Compute 6 as follows (reference 8) 


6 = a- tan 1 ^ i tan o J 


'(106) 


11. Compute qj as follows (definition of terms following equations (81)) 


_1 


$ = tan [m tan 6] 


(107) 


12. Compute 0 from equation (80), i.e. 
cr 


6 = tan 

cr ' 


■‘V 


(y + 1) (m-1 ) 
2m^ 


(108) 
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13. Test for type of solution; e.g. 

a . Tf _ I 6 I < 0 < |e I then "subsonic" solution 
1 cr 1 1 cr 1 

ap plies. Go to Step 14. 

h. If - 90 < 0 < )8 I or I 0 I < 6 < 90 then 

u - 1 cr 1 1 cr ' 

"supersonic" solution applies. Go to Step 29. 

14. Compute (3 as follows 



where 

0 = U/a 

U = — W . sin a . 
m A 


15. Compute £5^ as follows 



where 

V? = U/cos <J) . 

16. Compute C, D, E and F as follows 


C 


D 


v-i + 3 ~r 

Y+l Y + l 


m 


J tan 4> - £ (m-l) + 

1 sin (j> cos <f> 
Y + l J 


~ (m-1) [1 + (m-1) cos <£>] 


(109) 

( 110 ) 
( 111 ) 

( 112 ) 

(113) 


(114) 

(115) 


D 


(116) 
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17 . 


18. 


19 . 


E = 2 [ 1 " 77T m ] + 2(n ’“ 1) 


3 cos <j> 
w 


3 


F = ~~ [2 ( m - 1 ) sin <j> cos ij> ] 
8 


F = — F 
6 

w 


Compute a , b , c and d as follovis 


a = m 


CE + DF 

2 2 
C + D 


b = m 


CF - DE 

2 2 
C + D 


c = - D - F 
m 


d - > D* 
m 


Compute A and B as follows 


A = sec <f> + 2(m~l) cos $ + — (m-l) sin d> 

m 

b 2 

B = — (m-l ) sin $ 
m 


Compute $ , 6^ and <5^ as follows 


6 = - tan 1 


r _ 2 

U tan <j> 
8 2 


5 = t an 

s 


6 = tan 

P 


-1 


-1 


~B 

A 


c 8 - d tan d> 

w 

d 8 + c tan d> 

w Y 


(117) 

(118) 
(119 ) 

( 120 ) 

( 121 ) 

(122 ) 

(123) 

(124) 
(125 ) 

(126 ) 

(127) 

(128) 
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20 . 


21 . 


22 . 


22 . 


23. 


Compute S as follows 


* = ^ 


2 2 
A + B 


(129) 


Compute n (the subscript zero indicates that II is 
being evaluated at x = x = 0 ) as follows 


2 ~ 2 

n = ™ \/ c + d 
o m p 


(130) 


Compute the ratio of the refracted shear wave velocity 
amplitude to the incident shear wave velocity amplitude; 
i . e . 


r , ' 2 ~ 1 2 | -. 1/2 

U»s + w s I ] 

1 w . I 


= S sin o 


(131) 


Compute the ratio of the generated pressure wave velo- 
city amplitude at x = 0 to the incident shear wave 
velocity amplitude; i.e. 


~ 2 -' 2 . 1/2 


C 1 w + w I ] 
p p 


w , 


= H sin ct 
o 


(132) 


Compute the rms pressure just after the shock by us-ing 
equations (102) and (130); i.e. 

/2 y m | sec 


(6P) 


= £ 


where 


rms 


2 Y 


(y + 1) m - ( y 


4||n l! 

- 1 ) O 


(133) 


2 _ y-1 

Y+1 '“’A Y* 1 


(reference 8 or Equation 
(64)) 


(134) 


1 atm 


14.696 psia 


(135 ) 
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24 . 


25. 


2 6 o 


27 . 


28. 


Compute the sound pressure level just behind the shock 
i.e. at x = 0 , using equation (103) as follows 

( 6p ). 


SPL = 20 log 1 o — 
o u p 


riris 


(136 ) 


ref 


where 


-8 


ref 


= 0.29 x io psi. 


(137) 


Compute the surface wave number k for the pressure 
r sur r 

waves by using equation (96) as follows 


sur 


cos (a + 4> -5) cos 6 
cos (p 


k . 


(95) 


Compute the pressure wave frequency f from equation 
(99) as follows 


M A k a A 

f - 1 cos (a+0) 

2TT 


(13 8) 


Compute the ratio of surface pressure wave velocity to 
the local speed of sound from equation (101) as follows 
V p 

F sur _ cos <j» cos ( q+ 8 ) 


M 


A 


cos 0 cos ( a+<j> - 6 ) 


A a 


( 101 ) 


Compute the shock front displacement amplitude | <$x 


and phase 6 , , as follows 

^ shock 


6 x 


K cos 0 


.2 .2 
a + b 


(139 ) 


and 


{ sbock = tan ‘ 1 (a/£) - 


(140) 
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The steps that follow are for the "supersonic" solution 
and they correspond to the foregoing Steps (14) through (28). 

29. Compute 6 as follows 


1 - U 


30 


Compute as follows 


8 . 


w 




where 

W = D/cos <j) 0 

l i t r 

31. Compute C , D , E , F and G as follows 


C = 2 


Y-l 

Y+l 


m - 2 [1 + (m-l) cos <}0 


t 2 

D = (jn-1) [1 + (m-l) cos $ ] 


1 z f 3-Y 

E = (m-l) sm (J> cos $ - ^ 1 + — m tan <j> 


F -2 (m-l) sin <J> cos (j> 

I 

G = tan (p-<j)) 
where 


P = tan (1/8 ) . 

V7 


32. Compute a and c as follows 


a = m 


T t t 

C + G F 

i i t 

E + G D 


i i 

c = - D - F 
m 


(109 ) 

(141) 
(113 ) 

(142) 

(143) 

(144) 
(145 ) 

(146) 

(147) 

(148 ) 
(149) 
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33 . 


34. 


35 . 


36 o 


37 o 


38 . 


Compute A as follows 


A = sec 4> + 2 (m-1) cos q> + — (m-1) sin <t> 

m 


( 150 ) 


i 

Compute <£ as follows 


(j> = ] [ | <(> | - y 3 1 with sign of <j> 

or 

* = yjy I C U I - 4 3 | 

A A 

Compute S as follows 


A 

S = A /m 


(151) 


(152) 


Compute II as follows 


n 


c s m y 


m cos ( $-y ) 


(153) 


From equations (91) compute the ratio of the refracted 
shear wave velocity amplitude to the incident shear 
wave velocity amplitude; i.e„ 

2. 


cn: ♦ ;:*p 1/2 


= S sin cr . 


(154) 


From equations (91) compute the ratio of the generated 
pressure wave velocity amplitude to the incident shear 
wave velocity amplitude, i.e. 


~ 2 - t 2 i/o 

[ 1 w + w j 3 
1 1 


= JI 


sm a 


(155) 
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39 . 


40 . 


41. 


42 . 


43 . 


44. 


Compute the rms pressure by using equations (102) and 
(153); i . e . 


(dp) 

rms 


/2 7 m sec <p ** 

E TyTll Si - ( y - 1 ) ^ 


where p is given by equation (134). 


(156) 


Compute the sound pressure level by using equation (103) 
as follows 


SPL 


20 log 10 


(dp ) 


rms 


ref 


where 


p ^ = 0.29 x 10 
ref 


-8 


psi 


(15?) 


(137) 


Compute the surface wave number k 

sur 

waves by equation (96) as follows 


for the pressure 


k 

sur 


cos (a + <j> -6) cos 9 
i 

COS 4 


k 


(96 ) 


Compute the pressure wave frequency f from equation 
(99) as follows 


f 


*A k a A 

2ir 


cos 


( o+ 6 ) 


(138 ) 


Compute the ratio of surface pressure wave velocity to 


the local speed of sound from equation (101) as follows 


V 

p 

c sur 
a 


* 0 ^ 
cos d> cos (a + 0) r f A 

— — — — M . — — 

* As 

cos 0 cos (a+9 -6) 


( 101 ) 


Compute the shock wave displacement amplitude j<$x| as 


follows 
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e a 

k cos 6 


(158) 


A FORTRAN IV program to make these computations is given in 
Appendix C. 


Results 


The shear wave-shock front interaction will be illustrated for 
the same two Saturn V interstage flare angles, namely, 6 a 8° 58' and 
6 = 16° 44' . 

Appendix D contains computed data for the following cases: 

1.50 2.00 2.50 3.00 1.75 2.00 2.50 3.00 

A 

6 8.98 8.96 8.96 8.97 16.73 16.73 16.74 16.75 

The data column headings are defined as: 

THETA - , the shear wave incidence angle, m degrees 

PHI-PR -<{>', equation (126) or (151) 

WS/WA - shear wave amplitude ratio, equation (131) or (154) 

DEL S - 6 , equation (127) 

3 

WP/WA - pressure wave amplitude ratio, equation (132) or (155) 

DEL P - 6 s equation (128) 

(DP)RMS - as defined by equation (133) or (156) 

KSUR/K - k /k, equation (96) 

SUP 

VSUR/A - V /a, equation (101) 

P sur 

ADX/EL - as defined by equation (139) or (158) 

DEL SH - 6 s ^ oc ^ , equation (140) 
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It should be noted that there is at least one shear-wave incidence 
angle for each case considered that produces shock front resonance. 



CHAPTER III 


EXPERIMENTAL INVESTIGATIONS 

A model experiment was designed to study the interaction 
of sound waves and an oblique shock front. From the analysis 
of this phenomenon presented in Chapter 2 it Is clear that the 
sound field after the shock will depend on the initial Mach 
number, the turning angle (oblique shock angle) and the angle 
of incidence of the plane sound waves. The experimental sys- 
tem to be described used a fixed Mach number, Mj ~ 1,6, fixed 
incidence angle, 6 = 90 - a, and a variable turning angle, 6. 
The sound pressure level difference across the shock front for 
these conditions is shown in Figure 16 . 

Instrumentation and Equipment 
Model Wind Tunnel 

The low speed supersonic wind tunnel shown schematically 
in Figure 17 was built for these tests. Atmospheric air 
entered the tunnel by flowing around the horn, down a 24-inch 
long calming section and then through a 16-to-l smooth con- 
traction. After passing through the tunnel the air entered a 
vacuum chamber. The area ratio at the wedge (model) was 1.531 
but the wall boundary layers reduced this to an effective area 
ratio of 1.250 with a corresponding Mach number of 1.6. 
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We dge 

The wedge was mounted on a circular holder with the wedge 
leading edge at the holder's axis of rotation (see Figure 17 ). 

A half-degree incremented protractor coupled with a pointer on 
the holder indicated the wedge angle-flow turning angle 6-to 
within a quarter of a degree. 

A static pressure tap 0 o 0135-inch in diameter was provided 

0. 2-inch from the leading edge. 

Static Pressure Probe 

A circular probe holder formed the tunnel wall below the 
wedge (see Figure 17 ) and its axis of rotation was also at 

the wedge leading edge. As shown in Figure 18 , the static 

pressure tap was located 0.265-inch off center and, therefore, 
it was possible to measure the static pressure along a 180- 
degree circular arc. A 0-50 psi Heise Bourdon gauge was used 
as a pressure transducer. 

Acoustic Pressure Probes 

Two acoustic pressure probes were used. The one shown in 
Figure 17 consisted of a one-inch long 0.5 mm (0.009-inch 

1. D.) Bruel and Kjaer Type DB 0240 probe mounted on a BSK Type 
4134 condenser microphone cartridge. A vacuum tight extension 
was built to couple the microphone cartridge and a B£K Type 
2615 cathode follower. A BSK Type 2107 frequency analyzer was 
used to determine the root-mean-square sound pressure level. 
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The probe tip mated with a 0„009-inch hole cast into the 
probe holder thereby producing a smooth surface at the sensing 
port. The port, as shown in Figure 18 , was located 0.345-inch 
off the center of rotation, and 20.5 degrees behind the static 
pressure tap„ The acoustic pressure, therefore, could be mea- 
sured along 180 degrees of a circular arc. This path inter- 
sected the oblique shock front at right angles . 

It should be noted that both the static pressure and the 
acoustic pressure are sensed at ports that are normal to the 
flow direction and beneath a viscous boundary layer. Since 
there is a subsonic region next to the wall, the increase in 
fime-mean pressure due to the shock is felt upstream of the 
shock fronts In an attempt to circumvent the boundary-layer 
effects a second probe was developed. 

As shown in Figure 19 this probe extended from the calm- 
ing chamber to the vacuum chamber. It was a BSK Type DB 0241 
(1 mm O.D.) probe that was modified by plugging the tip, dril- 
ling a 0.0135-inch hole in the side of the tip and adding a 
length of tubing to the other end of the plug. The probe could 
traverse a one-inch path in the flow direction and its axial 
position was set to within 0 o 001-inch with a micrometer. The 
boundary layer above this port was much thinner and the acoustic 
pressure so determined is believed to be much more indicative 
of that in the flow. 

The same microphone cartridge and cathode follower were 
used with this probe but a different extension was required. 
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Sound Generator 


The sound field was generated by an Altec 802D driver 
loudspeaker in conjunction with the conical horn - cutoff fre- 
quency of 450 Hz - shown in Figure 17 . A Hewlett-Packard 

Model 206A audio signal generator acted as the electrical 
power source, A Marantz Model 9 power amplifier was required , 
to increase the power output to levels up to 30 watts. The« 
power to the driver was monitored by a Fluke Model 102R VAW, 
meter and the frequency was measured with a Hewlett-Packard 
Model 5212A electronic counter. 

Regulated 60-Hz } 117-volt power was supplied for all 
electronic systems by a Sorensen Model 2501 A.C. line voltage 
regulator . 


Experimental Procedure 


The 
used for 
1 . 

2 . 

3. 

4 . 


following is a synopsis of the experimental procedure 
taking all of the data: 

The electronic systems were turned on and allowed to 
stabilize for at least eight hours. They were cali- 
brated just before making a data run,. 

The vacuum system, was started and run for at least 30 
minutes before taking data. 

The pressure gauge zero was checked against a high' 
vacuum . 

The desired turning angle was set by rotating the 


wedge holder 



5 0 


6 o 


7 o 


8 . 


9 . 


10 . 


11 . 


12 . 
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The acoustic pressure port was positioned upstream of 
the shock £ront 0 

Power was supplied to the driver at the desired fre- 
quency. (The system could be tuned to obtain higher 
SPL's by making slight changes in the horn's axial 
location and in the driving frequency.) 

The power was varied until the desired SPL was ob- 
tained. (Since the flow noise level was quite high 
the impressed sound pressure was "viewed" in an eight 
percent band about the driver frequency by operating 
the B6K 2107 in the frequency analysis mode.) 

The driver frequency, current, voltage and power were 
measured . 

The flow field was traversed by either rotating the 
probe holder or inserting the linear motion acoustic 
pressure probe. 

Sound pressure levels were recorded for each angular 
or axial position by time averaging the output from 
the driver frequency band with a B£K Type 2^17 random 
signal voltmeter. (Averaging times of up to 100 
seconds were required when the level was close to the 
random amplitude background noise.) 

Static pressure was recorded at each angular position 
of the probe holder. 

Finally, background SPL's in the same frequency band 
were recorded. (The field was probed without power 


to the driver . ) 



Results 


Fluctuating pressure data were taken for stationary oblique shock 
fronts irradiated with sound waves. Since two different probes were used 
to obtain these data typical results will now be presented for each probe. 

Tunnel Wall Tap Data 

Figure 20 shows a typical fluctuating pressure -distribution with and 
without an impressed sound field. Assuming that the background pressure 
field and the impressed sound pressure field are not correlated, the 
sound pressure distribution is easily obtained. It is shown in Figure 21 
together with the static pressure distribution. 

It should be noted that the static pressure starts its rise some 
30 degrees before the predicted inviscid shock front position. This 
rather steep adverse pressure gradient undoubtedly thickens the boundary 
layer before the shock which in turn gives rise to a possible separated 
flow below the primary shock and curved oblique shocks at the corners 
formed by the intersection of the primary shock front and the upper and 
lower tunnel walls . Thus the unsteady pressure component was sensed 
beneath a very complex flow region and a reasonable comparison with the 
predicted acoustic pressure level distribution is not to be expected. 

The dip m the sound pressure curve prior to the primary shock front 
position was characteristic of all data taken. As shown in Figure 16 
attenuation of the sound pressure by shock interaction was predicted and 
this premature dip may be due to , at least in part , an interaction of the 
sound waves and the corner shocks. 

The rise in the pressure level near the shock front is to be expected 
since the shock is oscillating relative to the pressure transducer with the 
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same frequency as the impressed sound waves. However the level was 
expected to decrease away from the shock and be approximately 9 dB less 
than the sound pressure level before the shock. As shown in Figure 21, 
this did not occur and it can only be said that there was no appreciable 
change in level. 

In an attempt to measure the aco,ustic pressure through a much thinner 
and hopefully less complex boundary layer the side tap needle probe was 
installed and the following data were obtained. 

Needle Probe Data 

Figure 22 is a composite of the pressure data obtained with the 
needle probe. The parameter is the power supplied to the acoustic driver. 
The premature dip is evident for the high power run but it is apparently 
absent for the low power run. The increase in sound pressure level as 
the shock front is approached is much more evident with this probe. 
Increases of 17.5 to 23.0 dB from the initial level were obtained. A 
characteristic not found with the wall tap probe was the sharp drop in 
level after the maximum, a subsequent rise again, and then a gradual 
decrease in sound pressure level. 

This second dip is felt to be due to the relative position of the 
oscillating shock and the acoustic pressure tap. If the peak-to-peak 
amplitude of the shock displacement is less than the pressure tap dia- 
meter then the apparent acoustic pressure should drop to a relative 
minimum when the shock's mean position is in the center of the tap. 

Finally, it should be noted that ~the sound pressure level down- 
stream of the shock did not go below its inixial level and for the low- 
powered and non-powered cases the minimum levels were at least five 



76 


decibles higher than their initial levels. 

Since a large increase in the pressure level at the shock front 
was found for the case of natural pressure fluctuations at a frequency 
of 1818 Hz, a frequency analysis was performed to determine the effect 
of the shock front interaction with these fluctuations at other fre- 
quencies. To perform this analysis the pressure, as detected by the 
uncalibrated probe, was tape recorded for two cases: (1) probe at the 

shock front and (2) probe in the same position but with no shock pre- 
sent (the model was removed). Figure 23 shows the results and, for 
clarity, the difference of the two spectra is shown in Figure 24. 



CHAPTER IV 


DISCUSSION OF RESULTS 

Since there are two different phenomena treated in this study the 
discussion is divided into a section on sound wave-shock front interaction 
and one on shear wave-shock front interaction. Finally, an appraisal of 
the status of each of the theoretical analyses is presented. 

Sound Wave-Shock Front Interaction 

Analytical Study 

The first part of this discussion is based on the results given in 

Appendix B. Figures 8 and 9 will be helpful for visualizing the orienta- 
# 

I 

tion of the incident sound wave — as given by 8 (THETA), the refracted 
sound wave — as given by <f>' (PHIP) , and the generated entropy-vorticy 
wave — as given by 8 (BETA). In looking at these data note in particular 
the behavior of the ratio of the amplitude of the acoustic pressures in 
Regions 1 and 2 (PP2/PP1). 

Although the analysis was for a plane sound wave it will be more 
instructive to think of the oblique shock front being simultaneously irra- 
diated with equal amplitude plane sound waves with orientations from 
6 = - 90° to 9 = 90°. Since the governing equation is linear these fields 
will superpose. Look at the sound field downstream of the shock front. 

For the case of = 1.50 and 5 = 8.98°, the sound pressure is amplified 
by the shock when -90° <_ 0 < -40°, attenuated when -35° < 0 < 70°, and 

amplified again when 75° < 8 < 90°. Contrast this with the case of - 3.00 
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and 5 = 16.75° where there is amplification for all orientations but 
-5° < 6 < 20° and where shock resonance occurs at 0 « 53°. 

Since shock resonance was predicted for all but one of the cases con- 
sidered, it is interesting to compare the incident sound wave angle which 
produces resonance, 0 , with the angle of the flow before 'the shock, 

9 _ = 90 - a. 
r 


6 M. 9 6 

± r : 



"1 


r 

8.98 

1.50 

35.6 

— 

8.96 

2.00 

51.8 

82 

8.97 

2.50 

59.1 

68 

8.98 

3.00 

63.5 

57 



1 

v f 

r 

16.74 

1.75 

32.5 

82 

16.73 

2.00 

42.2 

74 

16.74 

2.50 

51.1 

63 

16.75 

3.00 

55.9 

53 


For relatively low Mach numbers the sound ray is almost parallel to the 

* 

shock front and moving leftward; however, at higher Mach numbers the sound 
ray aligns with the flow. Thus plane sound waves traveling in the flow 
direction would be strongly amplified for the case of = 3.00 and 
5 = 16.75°. 


Although the change in acoustic pressure across the shock front is 
independent of frequency, the peak-to-peak shock front displacement, |2a|, 


is inversely proportional to frequency since 


2a | = | PPAIlPIi 



Experimental Study 

The experimental study highlighted the difficulty of interpreting 
fluctuating pressure data from a supersonic flow. It can be said with 
certainty that the impressed sound field interacted with the oblique shock 
and caused significantly amplified pressure fluctuations near the shock 
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front . However , the shock process in this miniature system was too com- 
plex to permit a comparison of predicted and measured sound pressure levels. 

A decrease of about nine decibels was predicted but no appreciable change 
was found except in the vicinity of the shock front. 

The most interesting data to come from this study was that which 

showed strong amplification of the natural pressure fluctuations in the 

flow at the shock front. Figure 24 shows an amplification of about 26 dB 

\ 

for frequencies from 20 to 2000 Hz and definitely indicates that the shock 
motion contains all the frequencies present in the shockless flow. 

This points up a clear need for the accurate determination of super- 
sonic-tunnel sound fields if vehicle surface pressure data are to be pro- 
perly interpreted. The analyst must be able to separate tunnel sound effects 
from vehicle self .noise effects . 

Shear Wave-Shock Front Interaction 

The results given in Appendix D show an erratic behavior for most of 
the computed properties. In particular, the transition from the "super- 
sonic" to the "subsonic" solution looks suspicious for <£' and <S . 

Since the results are in doubt it will only be noted that shock resonance 
is predicted for all eight cases and that the results for positive 6 are 
similar to those of Ribner. 

Status of the Analytical Studies 

The analysis of plane sound wave interaction with "oblique shock fronts 
has been independently performed by two of the' authors. Although the resulting 
forms of the solution differed slightly the predicted data agreed. Thus it 
is felt that this solution is reliable. 



The shear wave analysis is not ready for application. There were 
two functional differences between Ribner's development and the one pre 
sented herein that have not been resolved. Further, an independent 
check of this analysis has not been obtained . 



CHAPTER V 


CONCLUSIONS AND RECOMMENDATIONS 

From the plane sound wave-oblique shock front interaction study it 
is concluded that : 

1 . There exists , under suitable conditions , a sound wave incidence 
angle such that shock wave resonance will occur. 

2. The amplification of plane sound waves by the excited shock is 
independent of the frequency of the sound. 

3. The amplitude of the shock front displacement is inversely pro- 
portional to the sound frequency and directly proportional to the magnitude 
of the sound particle velocity. 

4. The analysis appears to be reliable enough for application. 

5. Oblique shock fronts may be excited by an impressed sound field 
and/or by a natural sound field. 

From the shear wave-oblique shock front interaction analysis it is 
concluded that: 

1. There appears to be a shear wave incidence angle such that shock 
front resonance is produced. 

2. The predicted results seem to be erratic in behavior and are, at 
this stage of development, suspect. 

3. The analysis is not reliable enough for application. 

As a logical extension of the sound wave-shock front study it is 
recommended that : 

1. The analysis should be applied to a multiple shock field of finite 
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extent with the objective being the prediction of the far sound field 
associated with specific sound wave irradiation, 

2. An experiment should be designed to confirm the prediction of 
shock resonance. 

To put the shear wave -shock front analysis on a firm base it is 
recommended that • 

1. The study should be repeated by another analyst. 

2. An experiment should be designed to confirm the predicted 
downstream velocity and pressure field. 



APPENDICES 



APPENDIX A 

FORTRAN IV PROGRAM FOR SOUND WAVE - SHOCK FRONT INTERACTION 
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303 

304 

no 


m 

314 

317 


321 

333 

333 

335 

343 

3 52 
357 
371 

373 

3 76 
101 

103 

116 

122 

123 

126 

127 

131 

132 
146 


146 

152 


152 

154 

156 

lfil 

162 

211 

214 


PK0SR4M *AIN( INPUT ♦OUTPUT * TAPES* INPUT ♦TAPE6»aUTPUT) 


THE EFFECT OF PLANE SOUND WAVFS ON A ST AT TONARY OBLIqUE SHOCK 
T 1 IS THE ABSOLUTE TEMPERATURE IN REGION l IN DEGREE*; RANKINE 
T 1 35l«. 7 

551 IS THE SPEED OF SOUND IN REGION l 
PLFa ( f 1.4*53. 35 H»32.1 74 )«Tl 
SSlsSnRT ( PLF ) 

THE WORKING FLUID is dry <\IR 
G is THE RATIO OF SPECIFIC HEMS 
6*1.4 

SHRatrtM.O)/ (G-l.O) 

PIs3 . 1415926536 

,$1 IS THE MACH NUMBER IN REGION \ 

SIGO IS THE OHLlQUE SHOCK ANGLE IN DEGREES 
DELTAD IS THf FLaRE (TURNING) ANGLE IN DEGREES 
REAO(5»101) Si .StGOfDELTAD 
FORMAT (3F1 0.3) 

SIGRsSIGD*Pl/l30 f 0 

A^Sl «S1#(SIN(SIGR)') #(SIN(SIGR) ) 

RR2i s (DENSITY 2) / (DENSITY 1) 

R»2(3 ( (G*1,0)*A)/(2.0+ (G-1.0)»A) 

PP21 s (PRESSURE 2 ) / ( PRESSURE 1) 
pP2l»l .0 + G*A* (1.0- ( l.O/RRHl ) ) 

DELTAR=S1GR-ATAN( (SIN(SIGR)/C0S (S1GR) )/RR?l) 
DFLTAO=DELTAR*180.0/PI 

T? IS THE ABSOLUTE TEMPERATURE IN REGION ? IN DEGREE^ RANKINE 
T?sTl»<PP21/BR21) 

552 is the Speed of sound in region 2 

SJFs( <T2*1 ,M *53.35) #32. 174 
SS2aSORT(SJF) 

S? IS THE MACH NUMBER IN REGION 2 

S?3 S1«(C0S(SIGR)/CUS(SIGR-0ELTAR) )*(SS1/SS2) 

PHIR TS THE ANGLE PHI IN RADIANS 
PHIH 2 PI/2,0+0ELTaR-SIGR 

THE FOLLOWING aRE CONVENIENT ABBREVIATIONS 
R 3 HR2 1 
AXaSl*SlN (SlsR) 
akx=ax*ax 

SXsSS?/SS1 

SXXaSX«SX 

WRITE (6» I(i8) S1*DEL-TA0,SIGD,S2 

FORMAT <1H1///20X,4HM1 = ,F6.2 »?X .7HOELTA s ,Fb. 2 , 2X , 7Hs IGMA a, 
F6.2»?X*4HM2 a f F5,?»2x.9HTHETACP s*9X» 

9HTHETACN =/) 

WRITE (6»l09) 

FORMAT (20 X, 90 UH-) /20X»5HTHETA*3X.4HBETA.5Xi4HPHIPi5x» 

7upp2/PPl,4Xt6HSPL2/l,3Xt7HU02/U0l,3X»SHUEVW2,4X» 
5 HSPRcU»4x»6HAT/UP1i3x»6HPPAMPL/20X*90(1H«) ) 

THD IS THETA, ANGLE OF INCIDENT SOUND WAVE IN OEGREEs 
T HD a- 05 . 0 
THOMHD + 5.0 

IF (tHD.GT .90,0) GO TO’ 100 

THR=THO»PI/laO.O 

BETAR=ATAN( (Sl*STN{ rHR)»COS(STGR)«COS(PHIH) )/{<SIN(P hIR))»( 1*0 
+Sl*SIN{SIGR)*OOS(THn) ) ) ) 

rftad=betar*iho.o/pi 

SUB-ROUTInE for FINDING PhI-PrIME 

TRAM= (SS2/SS) )»SIN( fHR)/(l .0+Sl tt SIN(SIGH*THH) ) 


C 


0 



230 

0 

OFLO=20,04PI/1«Q,0 

233 

0 

NaO 84 

234 

0 

IF(THri) 111 ,112,113 

236 

111 0 

PHlpHa2.0»BETAft 

240 

0 

OFLPRs-1 .0«*PI/180.0 

243 

0 

go re 114 

243 

112 0 

RFTARsO .0 

244 

0 

PHlPHaO.O 

245 

0 

GO TO 116 

245 

113 0 

PHIPHs2„0#BEtAR 

247 

0 

DELPRsl *0*Pl/lflO.D 

252 

0 

GO TO 114 

2-52 

114 0 

N=N*1 

254 

0 

TEST=nELP 

255 

0 

pHlpHsPHIPR+nEUPR 

2.57 

0 

GAJMaTH'AN«{1.0 + S2ftCOS (PHIR<'PHTPR) > /COS (PHI PR) 

271 

0 

PHIPP=AT4N(GAMR) 

e?3 

0 

OELP=ABS{ (PHTPrt-PHlPP) /PHIPR) 

277 

0 

IF(oELP,LT. 0.0001) 00 TO 115 

302 

0 

IF(N.GT»200) 60 TO 115 

305 

0 

IF (^EKP.LT.TpST) GO TO 1U 

307 

0 

PHlPHsPh[PR-?,0«nEl-PR . 

311 

0 

DELpRsjOELPR/p. 0 

312 

0 

DELP=nELP*l.f) 

314 

0 

GO TO 114 

334 

115 0 

PH1pOsPHIPR*1 80.0/PI 

C 


END PHI-PRIME SUh-ROUTINE 

316 ' 

0 

Ajs (2.QM/5HR-1 .0) *C0$ (THR) /Si A (SIQR) ♦ (G-1 . 0-R».(G - \ , 0 ) /SHR^ * 

- 

1 

6 1 

337 

0 

61 3- <?,0*R/ShR- 1,0> *C0S (SlGR) /SIN (SIGR) 

351 

0 

Cl a- (2.0 *R/ShR- 1.0> /SIN (SlGR) 

360 

0 

A2*SIN(THR) «SIN (PHIH) ♦SIN (SISRJ-ACOS (PHIR) 4A1/R 

375 

0 

823SIN(SIGR-)*SIN(PHIR) +SIN {SlGR ) «COS (PHIR) »81/R 

412 

0 

C?=SIN(SIGR) *C0S(PH1R)«C1/R 

421 

0 

A3sSlN(THR) *COS <PHIH) •COS (SlGR) <*COS (PHIR) Ml 

435 

0 

83*6 IM ( SIGR )*COS (PH 1R) -COS (SlGR) /S IN ( pH I N) -COS ( SlGR ) *COS < PHI R> 


1 


457 

0 

Cls-COS(SIGR) *COS (PHIR) *Cl 

466 

0 

A4s( 2.0*G*AXX-G* (G-l .0) > /.( 2. 0#G«4XX- ( q- 1 . 0 ) ) + 


1 

4.0»G4&X«COS(THR)/(2,0#G»4xX-(G-1.0) ) 

507 . 

0 

B4s-4 .0«G#AXX tt COS (SIGR)/ (SIN ( S IGH ) * ( 2 . 0*GMAX- (,G-1„ 0 ) ) ) 

524 

n 

C 4s~4.0*G#AXX/<SlN<SlGfl)«{2.0»G»AXX-(G-l.n> ) ) 

536 

0 

A5sl ,0-2. 04 t (COS(THK)/4X)+0.54 (G-1.0) )*< (R/SHR)-1.01 

552 

0 

653-2, 0*( ( R/SHR ) -1 • 0 ) 

556 

0 

C5=-iS/SlN (SlGR) 

561 

0 

BXaSl^SlN (SiGR+THR) 

566 

0 

CXaSl^SlNTTHR) 

571 

0 

Al Is < 1 ,0+RX) »C2-CX* \2 

576 

0 

A 1 2s 1 . 0 

577 

0 

A 1 3 -jCOS (PhIR-PHIPR) 

603 

0 

A 14 = 0 • 0 

604 

0 

A^1=(1.0+6X)#C3-CX«H3 

611 

0 

A22=1.0*COS(PHlR-8ETAR)/SlN(PHlR«8ETAR> 

623 

0 

A?3s-SIN(PHIp-PHIPR) 

630 

0 

a24 = 0 .0 

6*3 1 

0 

A31=(l ,0+HX) *C4-CX tf rf4 

636 

0 

A 32 = 0 , 0 

636 

0 

A33-G»SSI«S1/SS2 

642 

0 

A 34 = 0 • 0 

642 

0 

A4ls(l .0+RX)«CS-CX*>i5 



85 


b47 

350 

b52 

f>53 

b54 

356 

357 
561 

377 

716 

736 

756 


342 

344 

346 

351 

353 


372 


3 74 


377 

LOO 

130 
L-30 

131 
L 32 


0 A42*0 # 0 

0 A43cSl/SX 

0 A44 a 1 . 0 

0 01 »a2 

0 0?=a3 

0 03*44*51 

0 D4»aS*51 

1 A3l*A?2#ai3-All*A32*A?3-4i3»A2i*Al2 

0 EST4«*<Ol*A22»A33 + 02*A32*4l3+03#A23*Al2« 

1 01 *A23*AJ2-U2*4l2*A33-U3*A22«Ai3) /nE T 

0 FST4R*(Dl*423 tt A3'l*0 2*411*433* D3 # A13*A 2,1- 

1 01*421 *433-02*43 l*A13-03*A 11*423) /uET 

0 (JST ARs <0 1*43 ?* a2\ *02*4 12* a31 *03*411*422- 

1 . 01*A3l*422-02*ALl*A32-r)3 <t A21*A12} /OET 

0 HSTaRs <-Ul*U2 1 *A32#A43 + A31 *Aa2*A23 + A41*a33*A22- 

1 A41*A32<»A23-A3HA22*A43-A21*A33*A42) + 

2 02* ( A 1 i*A'32*A43*Al?«A33«A41*A13#A42*A31- 

3 A4i*A32*Al3-A31«Al?*A43-An*A33»A4 2)- 

4 03* { All*A22 tt A4 3*Al2*A23*A41*A3.1*A42*A21« 

5 . A41*A22*A13-a21*A1?*AA3-A11*A23»A42)+U4*0ET) /DET 

C PPPVl IS THE HaTiO Of THE PERTURBED PRESSURE IN REGION 2 TO 

C THAT IN region 1 evaluated AT THE SHOCK FRONT 

0 PPP21 =USTAR»PP21/,SX 

0 RPPp2aABS (PPP21) 

C S D L? I = SPL2 - SPU 

0 SPL2la20 .n*ALOGlo ^RPPP2) 

C U 02 0 1 IS THE RATIO OF THE SOUmO particle VELOCITY IN 

C REGICN2 TO THAT IN REGION 1 EVALUATED AT THE SHOCK F^QNT _ - 

0 U02OUUSTAR 

C UEVW? IS THE ratio of the rms of the e-v velocity to THE 

C PRODUCT OF The UnPEKTUHQEo VELOCITY IN REGION 2* I.E. W2» A„ u . 

C THE ACOUSTIC MACH NUMBER IN REGION 1 

0 uEV’J2s (SIN (PmIR) *FSTAR) / (SQRt (2»0) *S1*C0S (SlGR) » 

1 SIN(SIGR~BET4R)} 

C SPRCO IS THE RATIO OF the AMPLITUDE OF THE ENTROPY 

C PFflrURRATION IN RFGION 2 TO THE PRODUCT OF THE CONSTANT 

C PRESSURE SPECIFIC HEAT AND THE ACOUSTIC PARTICLE VELOCITY 

C IN REGION 1 

0 SPRCUa-HSTAR 

c a tup l is the Ratio of the shock front velocity to the 

c sound particle velocity In region i evaluate at the 

C SHOrK FRONT 

0 ATUPl=-ESTAR« (lon+BX) 

C PAM^l IS THE RATIO OF THE PRODUCT OF THE PEAK-TO-PEA|< 

C AMPLITUDE OF THE SHOCK FRONT MOTION AND THE CIRCULAR 

c frequency IN REGION 1 TO THE scunu particle velocity in 

C RFGTCN I 

0 p AMPLa2 .O^ESTAR 

0 WRITE (6*116) THD.RETAOtPHTPO,PPP21*SPL2I*U02Ql*UEVW2, 

1 SPRCU»ATUP1»PAMP| 

116 0 FORMAT{2nX,F5.1»FR.3»F9.3 f Fl0,4fF10.3«Fl‘1.4»4F9,3) 

0 GO TO 110 

117 0 N*N 

FNO 


SED COMPILER SPACE 
300 



APPENDIX B 

SOUND WAVE SHOCK FRONT DATA 



THETACN 


K1 = 1.50 DELTA = 8.98 SIGMA = 54.44 M2 a 1,16 THETACP a 


THETA BETA 

PH1P 

PP2/PP1 

SPL2/1 

U02/U0f 

UEVW2 

SPRCU 

AT /UP 1 

PPAMPL 

-90.0 -41,551 

-89,977 

1.5481 

3,796 

’ 1.0528 

-.010 

.111 

-.030 

.466 

-85. n — 38.5Q2 

-87 # 3 0 9 

1.5561 

3.841 

1.0583 

.011 

-.268 

-. 032 

.271 

-30*0 -j-35.763_ 

-80.525 

. -1*4896 

3.461 

. .1*0131. .. 

__ .004 - 

- -.794. 

-.009 

- .051 

-75.0 -33.049 

-74.130 

1.4286 

3.098 

,9716 

-.007 

-1,289 

. 03U 

-.128 

-70,0 -30.439 

-68.066 

1,3704 

2,737 

.9320 

-.020 

-1.759 

.083 

-.277 

*65 .-0 — 2-T-.922 — 

.—62,27-4 

--■-1 »-3-l 29 

2*364 

-.8929 

— -—.034 

- — 2.207 

.146 

•-.403 — 

-60.0 -25,488 

-56.726 

1.2553 

1.975 

,8537 

-.047 

— 2 » 635 

.210 

-.510 

-55.0 -23.127 

-51.380 

1.1970 

1,562 

.8141 

-.059 

-3.041 

' .29? 

-.604 

. - - -50.0 -20.832 _ 

-46.211 

1.1378 

- 1,121 . 

. -.7738 - 

_ . 069_ 

-3.425 

,382 

-.685, 

-45.0 -10.595 

-41.189 

1,0776 

.649 

.7329 

-.076 

-3.785 

.472 

-.757 

-40.0 -16.408 

-36.308 

1.0167 

• 1 44 

.6915 

— • 080 

-4.121 

.564 

-.821 

-35. a —-14.265 

-31.537 

.9554 

-*397 

. . ,64-9?. . 

.-081- 

_ -4,430 

.650 

-.878 

-30.0 -12.159 

-26,862 

.8939 

-.974 

.6079 

- • 078 

-4.711 

.752 

-.928 

-25.0 -10,086 

-22.269 

.8327 

-1,590 

.5663 

-.073 

-4.961 

,846 

-.973 

. . ..-20 _ .-8.038- 

- 1 7-. 7 4 0 

- .7-723- 

-2,244 

*525?— 

— - .064 — 

— - -5. 1 80 - 

.937 

-1,014_ 

-15.0 -6.011 

-13.260 

.7133 

-2,935 

,4851 

— * 052 

-5.366 

1.026 

- 1 . 050 

-10.0 -3.999 

-8.820 

.6562 

-3,659 

,4463 

-.037 

-5.518 

1.108 

-1. 081 

*c n .-1.9 9J- 

-4.403- 

- *6017. 

- — 4-. 4 J-2. 

--*.41)92—. 

-.*020--- 

—5 .634 

- 1*186 

- - 1 . 109 — 

0.0 0.000 

0,000 

.5505 

-5.185 

.3744 

,000 

-5.713 

1.257 

-1.133 

5.0 1.997 

4.403 

.5032 

-5.965 

.3422 

,021 

-5.755 

1.321 

-1.103 

inn 9 . Qqg 


,4607.. 

-6 ,-7-32 

,3133 _ 

*.0.44 . 

-5.-758 - — 

--1*375 - 

. _ -1..169— 

15.0 6.0U 

13.260 

.4238 

-7.457 

.2802 

.067 

-5.722 

1.419 

-1.180 

20.0 8.038. 

17.740 

.3934 

-8,104 

.2675 

.091 

-5,647 

1.452 

-1.187 

, n in,ns6 

--22,269 

,-37 05. 

_*8*625 

*2520 

.114 - 

..-5,532 

1,472 

-U1B9„. 

30 1 0 12.1=9 

26.866 

• 3563 

-8,963 

,2423 

,.135 

-5.377 

1.477 

-1,18b 

35.0 14.265 

31.537 

.3523 

-9. 61 

.2396 

.154 

—5. 180 

1.468 

-1.174 

Art ft 16-4A8 -56.308 

. * 36 U 2 - 

-ft.nha 

.*2450 . 

*1-70. 

-4.941 

-.1.441 

.. -1,15S_ 

45.0 18.595 

41.109 

^ 3824 

-8.351 

,2600 

.181 

-4,659 

1.394 

-1.124 

50,0 20.032 

46.211 

.4218 

-7.498 

,2868 

.186 

-4.330 

1,324 

-1 , 000 

__ ssj 

51.380 

.402.9- 

--6.323 

,3284 — 

.184— 

— -3.951 - 

- 1.220 

-1*017. 

60.0 25.4=8 

56.726 

i5723 

-4.848 

.3892 

.173 

-3.513 

1.097 

-.928 

65.0 27.922 

62,282 

,7008 

—3 , 088 

.4766 

.150 

-3.003 

.922 

— .800 

7 n n 3n-43<3 


-.-0 fi.6 7_ 

. -1 , r)44 

,-603-1 

*112 

. -2.396 . 

. 685— 

- . 6 13 — 

75.0 33.049 

74.098 

W 1*649 

1.325 

.7922 

.056 

-1.642 

.351 

-.326 

80.0 35,763 

80.525 

1.6160 

4.169 

1,0990 

-.025 

-.622 

-.159 

.153 

85*3 38*592-. 

..87,184 

2.4376 

- 7.739 

1.657ft 

-.134 . 

.960 

-1.041 

. -1.054. 

90.0 41,551 

89.852 

3.3588 

10.524 

2.2843 

-.162 

2.634 

-2.104 

2.247 



Hi B 

2.00 DELTA * 6,96 

SlfiMfl a 

38.20 M2 

s 1.68 

THETACP a 

THETACN ts 


theta 

beta 

PH 1 P 

PP2/PP1 

SPL2/1 

U02/UO 1 

UEV92 

SPRCU 

AT/UP1 

PPAMPL 

-90.0 

-41,343 

-89.999 

1.4736 

3.367 

• 977 £ 

-.045 

-.500 

.101 

.354 

-85. o 

-38.352 

-87.454 

1.5280 

3 * 

1.0131 

-.011 

-.738 

.048 

.210 

- - _ -so . 0 

- -35, 50-0 — 

-80.594 

1.5159 

3.614 

1.0050 

-.000 

-1.126 

.007 

. ...040 _ 

-75 . n 

-32.773 

-74.077 

1.4989 

3.515 

.9937 

.004 

-1.531 

-.010 

- » 105 

-70.0 

-30.157 

-67.935 

1,4768 

3.386 

.9791 

.003 

-1.949 

-.00b 

-.229 

— 6S- r 0 

_ .-2-7 , 640 

6 2 « 0 91 - 

1.4486 

- -3,219. — 

.. ...9604 

— _ -.001 

-2.380 - 

. -.017 _ 

js ,338__ 

-60.0 

-25.2] 1 

-56. 5r>5 

1.4139 

3.008 

.9374 

-.008 

-2.818 

.056 

- * 435 

-55.0 

-22,861 

-51.140 

1.3726 

2.751 

.9100 

-.017 

-3.263 

.110 

-.521 

„ -50„ 0 

... —20 «580 . 

.-45.961 

1.3245 _ 

__ ,2.441 

-.8781 

-- 026- 

_ -3*709 _ 

. .177 

. -.599 

-45.0 

-18.360 

-40.947 

1.2700 

2.076 

.8420 

-.034 

-4.154 

.256 

— . 67 1 

-40.0 

-1 6. 1<?3 

-36.074 

l.?092 

1.6S0 

.8017 

-.042 

-4,594 

.345 

-.736 

— 35-. 0 . 

— -14., 073 . 

-3U321 ~ 

- - 1 , 1 426- . 

._ -1-156— 

-.757-5- 

— o-048 

- -5. 026 

,443 _ 

--.797- _ 

-30.0 

-11,992 

-26.667 

1.0706 

,592 

,7098 

-.052 

-5.447 

.549 

-.854 

-25,0 

-9,944 

-22.090 

.9937 

-.055 

.6588 

- . 0 54 

-5.852 

.661 

-.907 

. . ... .-20-.Q 

- *7-9 23™ 

.-17.600- 

--*9124 ... 

«-,J.96 

.6049- 

052 

- -6.238 

.777 


-15.0 

-5.924 

-13.154 

.8276 

- 1 . 644 

,5487 

-.046 

-6.603 

.898 

-1,004 

-10.0 

-3,940 

-8.740 

.7398 

-2.617 

,4905 

- , 035 

-6.942 

1.020 

-1,049 

s5.il. 

.-1,968. . 

... -4.J67 _ 

,6499 

_.-3.743 

— *-4 3 O.Q . 

Q2Q- 

-7., 254. 

- .1,144.. 

_«1.D92_ 

0,0 

o.ooo 

0.000 

.5585 

-5.n59 

.3703 

.000 

-7.535 

1.267 

-1.133 

5.0 

1.968 

4.367 

. 4666 

-6.621 

.3094 

.026 

-7.782 

1,388 

-1.171 

1 o , n 

8, 94D 

8., 740 

.3751 

-8,518 

.2487 

, 06H 

-7..S94 . 

1 ^ 0 ^ 

-1,209 _ 

15 ! 0 

5.924 

13.152 

^2848 

-10.909 

.1888 

.096 

-8.167 

U618 

-1.244 

20,0 

7.923 

17,600 

.1968 

-14,120 

.1305 

.141 

-8.299 

1.725 

-1,278 

5>8,0 

... . 9.. 9 A 4 

_22»09&. . 

- 4 131 

—19, np 6 

,0743 

, 1 9"? 


__ 1.825 

-1*310 

3 0 j 0 

11,992 

26,671 

9 1 — — * — 

.0320 

-29^900 

.0212 

.253 

-8.432 

K915 

-1 .340 

35.0 

14,073 

31.317 

— . 04 22 

-27.486 

— , 0200 

.323 

-8,429 

1.994 

-1.368 

.. 4 0^.0 

L6.193. 

_J36.J)74_ . 

. -_,1 089 .. 

-1 9.757 

.-,.0722- 

. .408 

-8.374 

2 . 06 Q__ 

-3.393 

45.0 

18.360 

40.955 

- , 1660 

-15.600 

-.1100 

.494 

-8.264 

2.111 

-1.414 

50.0 

20.580 

45.973 

-.2104 

-13.541 

-.1395 

,600 

-8.093 

2.143 

-1.429 

.. -55.0. 

_ .22.861 _ 

51,129 

7_4 — 

. -12*49.2 

_ — « ,15 7 4 . 

. ,.724 

_ ..—7 .850 _ 

.. 2.150 

..-1.435 . 

60.0 

25.211 

56,485 

-.2387 

-12.441 

-.1583 

.873 

-7.513 

2.123 

-1.425 

65.0 

27,640 

62.093 

-.1973 

-14.099 

-.1308 

1.053 

-7.038 

2.044 

-1.387 

JTIU-Q 

20-15 7. _ 

67,939 

— W ..46 73 — 

„.-23^ 473 .. 

-.0444 

. *4. 264 

__ -6.300 

1.859 

-1.282, . 

75.0 

32,773 

74.046 

.1276 

-9,693 

.2172 

1.367 

-4.814 

1.371 

— , 9b6 

80.0 

35.500 

80.500 

2.9763 

9.474 

1.9732 

-3,268 

2.941 

- 1 . 66b 

1,206 

J35.0- 

3B.35H _ 

..87,704 

-2,9341 

9.350 

-1,9453 

-169,830 

.-13.052 

4,942.. 

-3.697. 

<?0.0 

41.343 

89,93? 

-2,2993 

7,212 

-1.5244 

-6.139 

-10.917 

4.251 

-3. 306 



Ml 


THETACN a 


s 2,50 DELTA = 8,97 SlftMA a 30.89 M2 = 2,13 THETACP = 


THETA 

BETA 

PHiP 

PP2/PP1 

SPL2/1 

U02/U01 

UEVW2 

SPPCU 

AT/UP1 

PPAMPl 

-90.0 

-40.802 

-89.978 

1.5191 

3.632 

,9405 

-.049 

-1.029 

,160 

.293 

-85.0 

-37.720 

-88,439 

1.6116 

4.145 

.9978 

-.015 

-1.143 

.095 

.184 

-80. 0__ 

-34.8o7_ 

-80,738 

1.6195 

4.180 

- 1.0026- 

-.001 

. -1.436 

.014 

. .032 

-75.0 

-32.043 

-73.930 

1.6250 

4.217 

1.0061 

.OOfl 

-1.765 

-.035 

-.094 

-70.0 

-29.411 

-67.588 

1.6228 

4,205 

1.0047 

.012 

-2.131 

-.059 

-.204 

— 65 . A~~ 
-60.0 

-26.896 - 
-24.4P3 

6l-»6i)4 

-55.923 

- 1,6115- 
1.5905 

4 , 1 45 — 
4,031 

, 99-7 7- 

,9847 

,012 

.008 

- --2.532- 
-2.963 

- -».q61- 
-.042 

_ . 302 

-.390 

-85.0 

-22.161 

-50.504 

1.5595 

3,860 

,9656 

.003 

-3,421 

-.005 

-.470 

— 50.0 ■ 

-19.938 

-45.308 

1,5184 

. -3.628 

- .9401 

- . » 004 ■ 

. - -3.902 

.049 

-.544 

- 45 . 0 

-17.745 

-40.302 

1.4671 

3.329 , 

.9083 

-.012 

-4.401 

• 1 20 

-.613 

-40 ,0 

-15.632 

-35.454 

1.4056 

2,957 

,8702 

-.020 

-4,915 

.205 

-.678 

-35.0 

-13.570 

-30.746 

- 1,3341 

- 2.-504 - 

.825 q 

— -.028 

-5.441 

.304 

-.740 

-30.0 

-11.553 

-26.153 

1.2529 

1.948 

,7757 

-.034 

-5*974 

.415 

-.799 

-25.0 

-9.573 

-21-654 

1.1622 

1.306 

,7195 

-.038 

-6.511 

.538 

-.857 

-20, n 

-7.623 - 

-17.234 

1.0625 - 

-- -,152-7-—. 

,657ft. 

-.039 

__ -7.048 

.672 

- - -.913 - . _ 

-15.0 

-5.697 

-12.874 

, 0540 

-.409 

,5906 

-.037 

-7.583 

. ft 1 5 

-.968 

-10,0 

-3.708 

-8.558 

.8372 

-1.543 

.5183 

-.030 

-8,113 

.967 

-1.022 

_ mi 5-®-0 — 
0,0 

- -1.891 

o.ooo 

- -4.272__ 
U.000 

-- .7124 _ 
.5799 

— -2.945 
-4,733 

.44U- 

,359o 

. .-,01ft 
.000 

-0,634 
-9 . 145 

- 1,126 
1.293 

-1,077. 

-1.132 

5.0 

1.091 

4.272 

.4399 

-7.133 

.2723 

.026 

—9 , 645 

1.465 

-1.189 

1 0 * 

3 T 7 Oft 

. .8,-559 

,3924 

— 1 0 r h ft 1 

r 1 ft 1 o 

.061 

-10.133 

1_.64.4_. 

-1 .247 

15 !o 

5,697 

12.874 

.1371 

-17.261 

,0849 

.108 

-10.611 

1.029 

-1.309 

20.0 

7.6?3 

17.234 

-.0267 

-31.467 

-.0165 

. 169 

-11.081 

2,020 

-1.374 

, 0 

9,57.3 

2L..6S4 

- 1 ?Q G4 

-13,961 

-«JL24l 

.247 

-11.549 

2.219 

-1,446 . 

30.0 

11.553 

26,153 

-.3867 

-8.253 

-.2394 

«35n 

-12.025 

2.429 

-1.526 

35.0 

13.570 

30,750 

-.5902 

-4.580 

-,3654 

.487 

-12.527 

2.654 

-1.617 

_ 40. n. 

_ _15, 632 

35.451- 

. --.3194- 

-1.730 

-_,507-3. 

.673 

_ -13.087 

2,904 

-1.727 

45.0 

17.745 

40.302 

-1.0907 

.754 

-,675 3 

.942 

-13.765 

3.197 

-1 .867 

50.0 

19.910 

45,305 

-1,4378 

3.154 

-.8902 

1.359 

-14.683 

3.571 

-2.059 

-S5.JL- 

— 22,1 6 1 — 

. -50,51-0 _ 

_ -1,9433. 

5,7.71 

_si_203l 

2,098 

: -16.136 

4.117 

-2.357 _ 

60.0 

24.4g3 

55.904 

-2.8613 

9.131 

-1,7715 

3,741 

-19.032 

5,119 

-2.926 

65.0 

26,896 

61.604 

-5.6084 

14.97? 

-3,472? 

9.958 

-28.412 

8.157 

-4.679 

. 7.0. . 0 

-29 ..4 1 1 — 

_..67.572. _ 

. 4J,;>3 58 — . 

. -32. 7-] 7 - 

_ 26.7676 

.-157.509. 

, „145.780 

-46,360. 

26.836.__ . 

75.0 

32.043 

73.961 

3.3314 

10.453 

2.0625 

8,336 

4.235 

-1.878 

1.103 

80,0 

34 ,8n7 

80,613 

1.3357 

2.514 

.827n 

-.391 

-2.419 

.308 

-l 65 «, 

-«$73 oo 

85.0 

37.720 

,87.439 .. 

. .5923 

. -4,549 

.3667 

-.803 

___-4,607 

... 1.093 

90.0 

40,002 

89.853 

.2067 

-13.610 

.1292 

-.726 

-5.499 

1.46b 

-.932 



Ml a 3. no DELTA b 8,98 SIGMA n 26,47 M2 = 2,56 THETACP a 


THETACN ® 


THETA 

BETA 

PHtP 

PP2/PP1 

SPL2/1 

U02/U01 

UEVW2 

SPRCU 

AT/UP1 

PPAMPL 

-90.0 

-40.235 

-89.969 

1,5989 

4.076 

.9185 

-.047 

-1.527 

.209 

.248 

-85.0 

-37.050 

-89.099 

1.7194 

4, 7o7 

.9874 

-.016 

-1 ,545 

.127 

.162 

-ao..o_. 

.-04*068.. 

80* 08 7 

1*7459 

4,840.- 

- 4*0026- 

—-*000 

-1.733 

.018 

*026 .. 

-75.0 

-31,265 

-73.779 

1.7693 

4.956 

1,0160 

.010 

-1.980 

-.054 

-.087 

— 70.0 

-28,616 

— 67 » 2 o 1 

1,7829 

5,n23 

1,0238 

.016 

-2.286 

-.099 

-.186 

-.65-^0 — 

—2-6-.-1-0-3--- 

— — -61,066 

- -1.7855- 

5*035 _ 

— -L, 0254 ,010 

--2*645- 

-*-.-i 19 

«. 27 5- 

-60.0 

-23.709 

-55.293 

1.7762 

4,990 

1.0200 

.017 

-3*054 

-.117 

-.356 

-55.0 

-21.418 

-49,817 

1.7543 

4.882 

1.0074 

.013 

-3.508 

-.093 

-.430 

-50.0 

-19*21-7- 

-44. 6.03 

-1.7194 

. . -4*7jia _ 

,9874. 

.338- 

--4,004 

. _ 049 

— -.500 . 

-45.0 

-17,094 

-39.604 

1,6711 

4,460 

.9597 

.001 

-4.537 

.013 

-.566 

-40.0 

-15.039 

-34.791 

1.6093 

4,133 

.9241 

-.006 

-5.106 

.094 

-.630 

—35.0 - 

-4 3*041 

-30.133 

-- 1*5336 

3* 7 1_4 

-*3307 

— *0.1 4_ 

— -5,705 

*192 

. -*692 

-30.0 

-11.092 

-25.602 

1.4440 

3.192 

.8293 

-.021 

-6,334 

.307 

-.752 

-25.0 

-9,184 

-21.180 

1,3405 

2 , 545 

.7698 

-.026 

-6.989 

.437 

-.812 

-20.0 

-7.-308 

-16*845 

1.2227 

1.746 

. -* 7 02 1 . 

--- -*029 

- -7.668 

*584 

-.873 

-15.0 

-5.4*9 

-12.576 

1.0906 

.753 

.6263 

— . 029 

-8.371 

, 74b 

-.934 

-10.0 

-3.629 

-a. 357 

.9436 

- ,504 

.5419 

- * 026 

-9.096 

.923 

-.997 

—5 »j0 — 

--i.au. 

.-4.1.70 

. *7814 

—2* 1 43 — 

, .4.4 6.7-. 

-.. 316 _ 

-9*846 

1.115 

-1,0 63_. 

0.0 

0.000 

0,000 

.6027 

-4.398 

.3461 

,000 

-10.621 

1.323 

-1.132 

5.0 

1.811 

4.170 

.4062 

-7 , H?5 

.2333 

.025 

-11.430 

1.549 

-1.207 

10.0-- 

.. 3.629 

_ -8*.35-7 . 

...189.3 

- 14 .457 

.1 007 

, Q62 

_ . 1 2 * 280- 

- - 1*794 


15’0 

5 ,459 

12.576 


-25,694 

— . 029g 

.115 

-13.188 

2.063 

-1.382 

^0.0 

7.308 

16.843 

-.3238 

-9.794 

-.I860 

.190 

-14.180 

2.362 

-1.488 

25.0 

9.1P4 

21*180 

-.6370 

- .-3,911. 

- —*3658- 

.296 

-15.300 

. 2.702 

-1.615- 

30.0 

11.092 

25.606 

-1.0100 

.086 

— , 58 0 0 

.452 

— 16 , 625 

3.103 

-1.773 

35.0 

13.041 

30.129 

-1.4769 

3.387 

-.848? 

.689 

-18.295 

3.602 

-1 .981 

.,40.0 _ 

. _UL« 039 . 

..34.797 

.. -2.1084- 

6.479 . 

. -.-1*2107- 

— 1,078 

-20.599 

. 4.272 

-2.278 .. 

45.0 

17,094 

39.595 

-3.0715 

9. ?47 

-1.7638 

1.797 

-24,220 

5.294 

-2.754 

50. Q 

19.21 7 

44.591 

-4,8871 

13,781 

-2.8065 

3.451 

-31.284 

7.224 

-3.689 

55.0 

21.418 

49*805 

--10,3353 

-_.20.286 

-5,9351 

- 9,616 

-53.125 

13.033 

-6.571 

60.0 

23.7o9 

55,293 

1 18.9562 

41 .508 

68.3111 

-180.013 

474,254 

-125,173 

62.676 

65.0 

26. 1 03 

61,082 

8.3402 

18.424 

4.7894 

-80.149 

23.808 

-6 . 96b 

3.484 

-70*0 

- -28*6 \ 6 — 

67*-232 . 

-- 4.0282- 

.. -42*0-02 — 

- 2*313^ . 

-..4-*aS2- 

— 6.671- 

-2,382 

-- 1.197— 

75.0 

31,265 

73.779 

2.4509 

7.787 

1.4074 

.691 

.701 

-.724 

.367 

BO . 0 

34.06B 

60,887 

1,6001 

4,083 

.9189 

-.083 

-2,295 

,155 

— .080 

-85*0- - 

-37-..050 - 

. 88,099 . 

-1.0637- 

— - - *b36. _ 

.6108- 

-.295- 

— -3,992 

-.691 

-*3 65 . 

90.0 

40.235 

89,969 

,7366 

-2,655 

.4230 

-.331 

-4.775 

.970 

-.526 


a 1.75 OELTa = 16.74 SIGMA a 57,50 M2 a 1,09 THETACP a THETACN - 


THETA 

BETA 

PH IP 

PP2/PP1 

SPL2/1 

U02/U01 

UEVW2 

SPRCU 

AT/UP1 

PPAMPL 

-90.0 

-39.028 

-89,994 

2.2632 

7.094 

1.0884 

-.057 

,064 

-.015 

.500 

-85.0 

-35.5P3 

-89.916 

2.3261 

7.333 

1.116& 

,008 

-.209 

-.032 

.334 

— 

—-32.473— 

--R1.523 . 

2, 1390 . .. 

6 *6f)4 

- 1.0287 

- _ .013 

-.888 

-.007 

- .041 

-75.0 

-29.534 

-73.521 

1.9885 

5.970 

.9563 

,005 

-1.474 

.041 

-.175 

-70.0 

-26.848 

-66.368 

1.8693 

5.434 

.8990 

-.007 

-1.981 

.104 

-.335 

0 

.74 . 344 G,9.ft44 

1 ,.7667 

. 4_. ^43 

—♦8496 

.. 020- - 

-2,437-.. 

- .178 

-.461- 

-- 6 0 1 0 

-21.995 

-53,841 

1.6736 

4^473 

^8049 

-.032 

-2.855 

.260 

-.564 

-55.0 

-19.778 

-48.252 

1.5856 

4,004 

.7625 

-.043 

-3.242 

.349 

-.649 

-50^0.. 

_-17* 674 

_-42.997 

. 1.4999 

. .3.522 _ 

-- .J213 

--.051 

—3 .600 

.443 

-.722 

-45 .0 

-15.667 

-38.030 

1,4154 

3,01? 

,6807 

-.058 

-3.933 

.541 

-.785 

-40'. 0 

-13.742 

-33.296 

1.3310 

2.483 

.6401 

-.062 

-4.240 

.641 

- . 840 

-3S.J3- 

_j=3J_*8p6_. 

--26.756 

_i*p4b2 . 

. 1.9 1 2 

,5993 — 

-..063 _ 

-4.521 

.743 

-.890 

-70.0 

-10.088 

-24,374 

1.1610 

1.297 

.5563 

— .061 

-4,777 

.845 

-.934 

-25.0 

-8.337 

-20.124 

1 .0^54 

« o32 

.5172 

-.057 

-5,006 

.94b 

-.975 

-?0 0 

-6. 6js5- 

„-l5*981 - 

...989.6- 

.-.091 

♦4759 . 

-.081 

«*5#208 

- 1.045 

-1.012.. 

-15> 

-4.943 

-11.916 

^9039 

-.878 

.4347 

-.041 

-5.382 

1.141 

- 1 .046 

-10.0 

-3.2°^ 

— 7 .913 

.8187 

-1.738 

.3937 

-.030 

—5 • 528 

1.234 

-1.078 

-5-.-0 

— -.l-*Jb-3.8-- 

_ .-J..946 — 

7344 — 

— -2 * 6B 1 — 

, 353? 

-..01-6 — 

— -5 * 644 

1.322 

■ -1.107 . 


0,0 

0.000 

0.000 

.6516 

-3.721 

.3133 

.000 

-5.730 

1.405 

-1.135 

5.0 

1.638 

3.946 

.5708 

-4,670 

.2/45 

.018 

-5.7R7 

1.482 

-1.161 

in,n 

3,?G4 

7.913 

.4.926 

,-iu LSO_ 

.2369 

. 037™ _ 

_*.s,ai3- 

1.552 . 

-1,186 

15^0 

4.943 

11.916 

.4177 

-7,582 

,2009 

.057 

-5.809 

1,613 

-1,209 

20.0 

6.625 

15.901 

.3468 

-9.199 

.1660 

.078 

-5.774 

1.667 

-1.231 

?5,0 

8l.33.7_ 

20.124 

,2805 

— 1 1 . n 4.1 

*33 4 q 

. 1.0 o._ 

-5.710 

1*711 . 

-1,25.1. . . 

30.0 

10.088 

?4 . 374 

.2198 

-13.161 

.1057 

.122 

-5,615 

1.745 

-1.270 

35.0 

11.886 

26.756 

.1655 

-15,625 

.0796 

.143 

-5.491 

1,769 

-1.288 

40.0 

13.742 

33.296 

.1385 

-18,500 

- t -D 5 72 

.164 . 

-5,337 

.3*732 . 

-1.3 03 . 

45.0 

15.667 

38.030 

.0814 

-21.790 

.039i 

,183 

-5.154 

1.782 

-1.316 

50.0 

17.674 

43.005 

.0553 

-25.150 

.026ft 

.200 

-4,941 

1.768 

-1.324 

. . .n . . 1 . . 5 5 . 0 . 

— 19. 778. „ 

__4B.244„ 

„ *0440... 

..-2 7., 124 

,0212 - 

.215 . 

-4.696 

1.737 . 

-1.328 .... 

60.0 

21.995 

53.833 

.0541 

-25,332 

.0260 

.225 

-4,414 

1.68b 

-1.321 

65.0 

24.344 

59.844 

.0993 

— 20 , o57 

. 047r 

,228 

—4 ,081 

1.605 

-1,296 

70 . 0 

26.848 66.384 

.? 1J2 

-1 3.X57 

,3.0.45 

. 21-9. 

.-3.658 

1.466 .. 

.-1,228 

75.0 

29.534 

73.568 

.5635 

-4,982 

.2710 

,179 

-2.973 

1.160 

-1.013 

80.0 

32.433 

81.616 

3.6647 

11.281 

1.7624 

-.194 

.795 

-1.122 

1.028 

. JB5,0_ 

_35.5S.3_ 

B9 .156 

-1 .9669 — 

_ _.5,976 _ 

-♦9459 

. ,455... 

- -5*3 87 

. -2,849.. 

-2*759. .. 0 

90.0 

39.028 

89.932 

-1 , 5364 

3.730 

-,738q 

* 360 

-4.363 

2,512 

-2.590 



I 


THETACN * 


M 1 = H*00 DELTA a 16.73 SIGMA a 47.80 M2 = 1,37 TriETACP = 


theta 

BETA 

• M ^ jNgJtf pp Ml 

PH IP 

PP2/PP1 

SPL2/1 

U02/U01 

-90.0 

-38,986 

-89.988 

2.1540 

6. 665 

1,0289 

—85 . 0 

-35.530 

-89,936 

2.2530 

7.055 

1.0762 

_«n n 

-32.374 

_-8l*537 

?. j 435 

6..622 

1.0239- 

-75.0 

-29.471 

-73.512 

2.0467 

6,221 

.9776 

<•70*0 

-26.7*4 

-66,333 

1,9651 

5.867 

.9386 


-i>5 ^0 -2A-ZPA m$9 . ao 2 

-60.0 -21.933 -53.787 

-55.0 -19.719 -48.195 

-JL7, 6 1 9 -42.940 _ .. 

-45,0 -15.616 -37,974 

-40.0 -1 3.6g5 -33,242 

_ -35-jQ. -li-.844.-_. -28.704 .. 

-30,0 -10,052 -24.330 

-?5.0 -8.307 -20.086 

-20,0 -6.601 -IS, 94a 


-15,0 

-10.0 


15.0 

20.0 

-2S.il — 

30.0 

35.0 

-4 0.-0 

45.0 

50.0 

55.0 

60.0 
65.0 

70.0 

75.0 

80.0 

85 . 0 

90.0 


-4.925 -11.893 


-3.271 
.-JU-632__ 
o;ooo 
1.632 
3 .2 7 1 - 
4.925 
6.601 
_a.-30-7 — 
10.052' 
11.844 
-13-6-95— 
15.6)6 
17.619 
-19._7l-9 „ 

21.933 
24.280 
26.-7.e4 — 
29.471 
32,374 
-35.530. 
■38,986 


-7,896 

-3.,-932_ 

0.000 

3.938 

7 .896 

11.893 

15.948 

20.-086- 

24.330 

28,704 

33.242- 

37,966 

42.940 

48*183- 

53.803 

59,810 

66-.318- 

73.567 
6l.49e 
—89. 061. 
89.972 


1.8159 
1.7408 
— 1.6627 
1.5808 
1.4944 
-3. 4035 . 
1.3080 
1,?080 
. .4, V03B— 
.9959 
,8845 

.77 0-1 

.6533 

.5342 

— ^-..44-35 — 
.2912 
.1676 

.-04-25— 

- . 1)845 
-.?148 
.35-04— 
-.4957 
-.6590 
— -.-8577- _ 
-1.1348 
-1.6209 
--3.-01-54 — 
29,6579 
1.7268 
.6640— 
.3197 


5-.S28- 

5.1 82 
4,815 
. 4,4.16- 
3.977 
3,490 
. _2.944 . 
2.332 
1,641 
.858 

-. i36 
— 1 * 066 
_ -2. 269- 
-3,698 
-5.445 
--7. 67-1- 
-10.715 
-15.513 
■,27.-424- 
-21.459 
-13,361 
.9 , 1 r>9 
-6,096 
-3,622 
-* 1-,333_ 
1,098 
4,195 

9.537- 

29.443 
4,745 
3 .556 
-9,904 


. 3026 - 

,8674 

.8315 

..1942. 

.7551 

.7138 

.67.0.4. 

. 624r 
.5770 

.S2-I3- 

.4757 

.4225 

.3120 

,255 2 

^1 SlS- 

.1391 

.0801 

.02-03- 

-,0404 
- , 1026 
— -...1-614- 
— . 2368 
-.3148 
— -.4097 

-.5420 

-.7742 

--1-.4-403- 

14,1663 

.8248 

-..3172 

,1527 


UEVW2 SPRCU 

» — «m m urn «• «. m _ mm »■» • ••■ 

-.060 -.332 

—.006 — ,518 

.609 -1,047- 


.011 
• 006 
--. 001 - 


-1,544 
-2'. 0 04 
_-2.438 


AT/UP1 

5 — rus 
.068 
.026 

„ -.001 

.006 

.036 

— .083 — 


Ol 0 

-2.855 

• ] 45 

— .503 

-.019 

-3,257 

,220 

-.586 

- . 02S __ 

_ -3.644 

„.304_ 

— . 658. 

— . 036 

-4.017 

.397 

-.723 

-.042 

-4.375 

,497 

-.781 

-.046— 

— *4.716. 

.602-. 

„ -.835- 

-.047 

-5.040 

.713 
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- 2.221 

2.315 

2.417 

—2.543 

2.724 

3.064 

—4.110 
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THETA 

BETA 

PH I P 

PP2/PP1 

SPL2/1 

U02/U01 

UEVW2 

SPrtCU 

AT /UP 1 

PP4MPL 
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-30.0 **9,532 -23.676 1.6179 4.179 ,6983 -.027 -5.427 .570 -.821 
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-15.0 -4.659 -11.543 1.2006 1,588 .5182 -.027 -6,641 .983 -.976 

-10.0 -3.093 -7.661 l.n381 .325 ,448i -.021 -7.034 1.136 -1.029 

---5.0 -1.543 - -3.820 .8634 -1.276 .-37-27 -*-04-2 -7.420 - --1.296 - --1.082 

0.0 0.000 0.000 .6758 -3,404 ,2917 .000 -7,800 1.463 -1.139 

5.0 1.543 3.820 .4741 .6.483 .2046 *017 -6.176 1.636 -1.198 

-WU-0 3.0-93 7.661 .2566 *-U.*£.i-6 -.41-07 038 -8.551 1. 821 «-l-263 

15.0 4,659 11,545 ,0202 -33,876 ,0087 .066 -6,932 2.016 -1,335 

20.0 6.248 15.489 -,?395 -12,413 -,103a ,100 -9.327 2,224 -1,416 

.25.0 7.869 19,528 --,5304 -5. 5^8 -.2289 .143 -9.752 - 2.451- . -1.511. 

30.0 9.532 23,673 -.8652 -1,258 -,3734 .198 -10,232 2,709 -1,626 

35.0 11.245 27.966 -1,?666 2.053 -.5467 .268 -10.811 3,013 -1,771 

40_a 13,070 — 33,439 .-1*7780 -4,958 -.Z6X4 .364- -11.572 . 3.397-___ -1 .567 

■>5.0 14.871 37.133 -2,4918 7.930 -1.0755 .502 -12.694 3.933 -2.256 

50.0 16.8U 42.075 -3.6426 11.228 -1.5722 .729 -14,622 4.801 -2.744 

55-0. L8,857 47,340 _ -6.0503.: 15.836 -2*6114 1U98._„18.914 . .6,632 . -3.796 _ _ 

60.0 21.010 52,997 -15.8699 24,011 -6,8496 3,073 -37.251 14,163 -8,163 

65.0 23.351 59.140 24.7208 27.861 10.6697 -4.537 40.071 -17.103 9.9B2 

0- r 6 25*8 *9 65*586 6 r 539-7 —16.,4-U 2*8226 -1*057 5.900— .-3. 144__— 4.868— . _ 


75.0 

28.558 

73,366 

3,5309 

10,958 

1.5240 

-.401 

.527 

— , 862 

,S25 

80.0 

31.517 

82.034 

2,2433 

7, *18 

.9682 

.006 

-1.574 

, 095 

-.059 

-8S*Q — 

_34.7J76*. 

. 89.552 - 

1*581.1 

-3.979 

_ _ ,6824— 

— .542- 

-2.476. - 

- .553 

_ -.359 

90.0 

38,395 

89.914 

1,2606 

2.011 

.5441 

6.390 

-2.710 

. ,706 

-.479 
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THETA 

BETA 

PH1P 

PP2/PP1 

SPL2/1 

U02/U01 

UEVW2 

SPRCU 

AT/UP1 

PPAMPL 

-90.0 

-37.805 

-89.985 

2.5187 

a. 024 

.9665 

-.047 

-1.302 

.174 

.235 

-85. n 

-33,9a2 

-89.965 

2.6922 
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1.0331 

-.012 

-1.324 
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.160 

^»80-. 0 — 

--30.598— 

--82,. 695 

- 2,7065 - 
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- 1.0386 - 

-- .013 

-1.561 

* .009 .. 

- -.016 

-75.0 

-27.570 

-73.390 

2.6624 

8.505 

1.0217 

.023 

-1,901 

-.061 
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-70.0 

-24.836 
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1.0077 
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-2.253 
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FORTRAN II-D PROGRAM FOR SHEAR WAVE - SHOCK FRONT INTERACTION 
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C 

C 

c 


0 

0 


0 


100 0 
101 


c 

c 


0 

0 


c 

c 

c 


0 

0 


c 

c 


0 

0 

0 


102 

1 


103 


1 

2 


C 

0 

104 


C 


0 


105 

106 • 
C 

107 0 


THE EFFECT OF SHEAR WAVES UN A STATIONARY OBLIQUE SHOCK 
T 1 IS THE ABSOLUTE TEMPERATURE IN REGION 1 IN DEGREES RANKINE 
T 1 = 5 18 . 7 

SSI IS THE SPEED UF SOUND IN REG I UN 1 

551 = SQRTF (1.4* 53. 35*32*17 4* Tl) 

THE WORKING FLUID IS DRY AIR 

G IS THE RATIO OF SPECIFIC HEATS 
G= 1.4 

SHR= ( G~ 1 . )/<G+l. ) 

Pl=3. 1415926536 

51 IS THE MACH NUMBER IN REGION 1 

SIGD IS THE OBLIQUE SHUCK ANGLE IN DEGREES 
DELTAD IS THE FLAKE (TURNING) ANGLE IN DEGREES 
READ 101, SI, SIGD, DELTAD 
FUKMATI 3F10.3) 

S1GR=SIGD*PI/180. 

AX=S1*SINF(SIGR) 

AXX=AX*AX 

RR21= (DENSI TY 2>/<DENSITY 1) 

RR 21= ( (G+l. )*AXX)/(2. + (G~l.')*AXX) 

PP21= (PRESSURE 2 ) / ( PR ESS UR E 1) 

PP21=1.+G*AXX*( l.-( 1./RR21 ) > 

DEL T AR = S IGR-ATANF ( ( S 1 NF ( S I GR ) /COS < S I GR ) J/RR21 ) 
DELTAD=DELTAR*180./PI 

T 2 IS THE ABSOLUTE TEMPERATURE IN REGION 2 
T2=T1*(PP21/RR21) 

552 IS THE SPEED UF SOUND IN REGION 2 
SS2=SQRTF ( 1.4*53. 35*32. 174* T2) 

52 IS THE MACH NUMBER IN REGION 2 
S2=S1*(C0SF(SIGR) /COS F ( S I GR-DEL TAR ) )*( SS1/SS2 ) 

R= { D ENS I TY 2 ) / ( DENS I TY 1)=M 

R=RR21 

THECR IS THE CRITICAL VALUE OF THETA 
THECK=ATANF( SQRTF ( ( (G+l. )*( K-*l. ) )/( 2.*R*R ) ) ) 

THECD=THECR* 180 . /P I 

PRINT 102, SI, DELTAD, SIGD, S2,THECD 

FORMAT) 1H1, 1H0, 1H ,10X,4HMA = , F6. 2 , 2X , 7HDE LT A =,F6.2,2X, 

7HSIGMA = » F6. 2, 2X, 3HM = , F5. 2 , 2X , 10HTHET A CR =,F8.2/) 
PRINT 103 

FORMAT ( 10X,88( 1H-) / 1 1 X , 5HT HETA , 2X , 6HP HI -PR , 3X , 5HWS /W A, 2X, 

5HDEL S,3X,5HWP/WA,2X,5HDEL P ,4X , 7H ( DP ) RMS , 2X , 6HKSUR / K, 
2X,6HVSUR/A»2X, 6HADX/ EL,1X,6HDEL SH/10X,88( 1H-)) 
THDN=-THECD 
THDP=+THECD 

THD IS THETA, ANGLE OF INCIDENT SHEAR WAVES IN DEGREES 
THD=-9(J. 

THD= THD + 5 . 

T HK=THD*P 1/180. 

PHIR IS THE ANGLE PHI IN RADIANS 
PHIR=ATANF ( R*SI NF ( THR ) /CUS F ( THR ) ) 

T EST 1 = TH D-8 5 . 

TEST2= THD- THUN 
TEST3=THD-THDP 
IF ( TEST1 ) 105,105,100 
IF ( TEST2 ) 108,108,106 
I F ( TE S T3 ) 107,108,108 

THE FOLLOWING IS THE SUBSONIC SOLUTIUN 
UBL=( AX* SSI ) / (R*SS2) 
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BE TAL= SQR TF ( l.-UBL*UBL> 

WBL=U8L/C0SF (PHIR ) 

BE TWL= SQR TF { l.-WBL*WBL) 

CL=( 1.0*SHR + ( (3.-G)/(G+l.))*R)*(SINF(PHIR)/C0SF(PHlR))- 
1 ( { R — 1. )*(K-1. )+2.*<K-l. ) /(G+l. ))*$ INF (PHIR >*CUSF( PHIR-) 

DL = ( BETWL*( R-l. ) /( B ETA L* BETA L ) ) *1 1. + l R- 1 . > *CUSM PH I R ) * 

1 CUSF(PhlR)) 

DPL=(BETAL*BETAL/BETWL)*DL 

fcL=2.*( l.-R*SHR)+2.*(R-l. )*(BETWL*COSF(PHIR )/BbTAL)* 

1 (BETWL*C0SF(PH1R)/BETAL) 

FL=( BETWL/ ( BETAL*BETAL) > *( 2 . * ( R- 1 . ) *S I NF ( PH I K ) * C USE (PHIR ) ) 
FPL= (BETAL*BETAL/BETWL)*Fl 
AWL=R*(CL*EL+DL*Ft)/<CL*CL+l)L*DL) 

BWL=R*(CL*FL-DL*EL ) / l CL*CL+DL*DL ) 

CWL=AWL*DPL/R-FPL 

DWL=BWL-DPL/R 

AL = 1 . / COSF (PHIR) +2.* {R-l. ) *COSF ( PHI R ) +A WL* ( R- 1 . )=MR-1. )* 

1 S INF ( PH I R ) /R 

Bl=BWL*(K-l. )*(R-i. )*SINF{ PHIR) /R 
C PHIPL IS THE ANGLE PHI-PRIME IN RADIANS 

0 PH I PL=~AT ANF ( ( UBL*UBl*SI NF ( PH I R ) ) / ( BETAL^BETA L*COSF (PHIR) ) ) 

C DSRL IS THE SHEAR WAVE PHASE 

0 DSRL=ATANF(-BL/AL) 

C DPRL IS THE PRESSURE WAVE PHASE 

0 UPRL=ATANF ( ( CWL*B tTWL-DWL*S INF ( PHI R ) /COSF ( PHI R ))/( DWL*BETWL+ 

1 CWL*SINF(PHIR ) /CCJSF( PHIR) ) ) 

ST AK L= SOR TF (AL*AL+BL*BL)/R 
PISTL=SQRTF (CWL*CWL+DWL*DWL)/(R*BETAL) 

C WSRWL IS THE SHEAR WAVE AMPLITUDE RATIO 

0 WSKWL = STARL#SINF( SIGH ) 

C WPRWL IS THE PRESSURE WAVE AMPLITUDE RATIO 

0 WPRWL=PI STL*SINF( SI GR ) 

DPKML IS THE RATIU OF THE RMS OF THE PERTURBATIUN PRESSURE 
TU EPSILON - IN PS1 (PRESSURE IN REGION 1 ASSUMED TO BE 1 ATM, 
I .E. P1=1A.696 PSIA 

0 DPRML=SORTF ( 2. ) *G*k*ABSF ( PI STL ) *PP21*1A .696/ 

1 (ABSF( COSF (PHIR) ) * ( < G+ 1 . ) *K-( G- I . ) )) 

SSUKL IS THE RATIO OF THE SURFACE WAVE NUMBER' TU THAT OF THE 
INCIDENT SHEAR WAVE 

0 SSUKL^COSF ( SIGk+PHIPL-DELTAR )*CUSF( THR ) /COS F ( PH I PL ) 

VSUKL IS THE RATIU OF THE SURFACE PRESSURE WAVE VELOCITY TO 
THE LOCAL SPEED OF SOUND 

0 V SURL = COSE (PHIPL) *COSF( SI GR + THR) *S1*SS1/ (COSE (THR )* 

1 COSF (SIGR+PHIPL~DELTAR)*SS2) 

DXWPL IS THE RATIU OF THE SHOCK FRONT DISPLACEMENT AMPLITUDE 
TO EPSILON TIMES THE WAVELENGTH OF THE INCIDENT SHEAR WAVE 

0 DXWPL=SOKTF ( A WL* A WL+BWL*BWL ) / l 2. *P I*COS F ( THR ) ) 

C DSHRL IS THE SHOCK WAVE PHASE 

0 DSHRL-ATANF(AWL/BWL) 

KD= 180 ./PI 
PHIDL=PHI PL*RD 

DSDL -DSR L*RU 
DPDL=DPRL*RD 
DSHDL=DSHRL*RD 

PRINT 109, TH0,PH1DL, WSRWL, DSDL, WPRWL, DPDL » DPRML > SSURL , VSURL , 

1 DXWPL, USHDL 

109 FORMAT ( 10 X , F6 . 1 , FB . 2 , F«. 3 , F 7 . 2 , F 8 . 3 , F7 . 2 , E I 1 . 4 , 4F 8 . 3 ) 

GO TO 10 A 

THE FOLLOWING IS THE SUPERSUN1C SOLUTION 


C 
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108 0 UBG={AX*SS1)/(R*SS2) 

BE TAG=SQRTF ( l.-UBG#UBG) 

WBG=UBG/C0$F ( PHIR ) 

BE TWG = SGRTF ( WBG*WBG-1 . ) 

C PG=2 * *{<*> SHR-2 . * { 1. + ( R-l • ) *C0S F ( PH I R ) *C0S F ( PH I K ) ) 

UPG= t K— 1 . )*( 1. + <K-1 . )*COSF ( PHIR) OS Ft PHIR ) ) 

EPG=(R-1. )*(R-l. )*SINF( PHIR )*C0SF( PHIR)- ( 1.+ 

1 .( 3. -G ) *R/( 1. +G ) ) *SI NF ( PHI R)/CUSF( PHIR) 

F PG-2 . * (R-l. )#SINF(PHlR)*C0SF{ PHIR) 

SUR IS THE MACH ANGLE IN RADIANS 

0 SUR=ATANF ( l./BETWG) 

GPG-5 IN F ( SUR-PHIR) /COSFI SUR-PHIR ) 

AWG=R*(CPG+GPG*FPG> /(EPG+GPG*UPG) 

CWG= AWG*DPG /R-F PG 

AG = 1. /COSF {PHIK)+2.*(R-1. ) *CUS F ( PH I R ) + AWG* ( R- 1 . )*{R-1. >* 

S INF ( PH I R ) /R 

PHIPG=(PHIR/ABSF( PHIR) )*ABSF< ABSF{ PHIR)-SUR ) 
PHIDG=PH1PG*18U./PI 
ST AKG = AG/R 

PISTG=CWG*SINF ( S UK ) / ( R-CUS F l PH I R-SUR ) ) 

WSRWG=STARG*SINF| SIGR) 

DS0G=0 . 0 

WPRWG=PISTG*SINF( SIGR) 

DPDG=0 • 0 

L)PRMG = SQRTF ( 2. ) *G*R* ABS F ( P I STG ) *PP2 1*14 . 696/ 

1 ( ABSF (CQSFIPHIR) )*( (G+l. )*R-(G-1.) ) ) 

SSURG = COSF { SIGR+PHI'P.G-DELTAK)*CUSF< THR ) /COS F ( P H I PG ) 

VSURG=COSF (PHIPG )*COSF( S I GR+THR ) *S 1* SSI / 

1 ( COSF(THKKCOSF( S I GR+PH I PG-UELT AR ) *SS2 ) 

OXWPG^ABSF ( AWG) /( 2 . *P I *CUS F < THR ) ) 

DSHDG=~90.*( THR/ABSFfTHR) ) 

PRINT 109, THD,PHIDG, WSR WG, DSDG, WPRWG ,DPDG ,DPRMG , S SURG , V SURG , 

L DXWPG , DSHDG 

GO TO 104 
110 N=1 

CALL tXIT 
END 

00009 5 1870000 -3 
00019 14000000-1 
00029 5335UOOO-2 
00039 32 1 74000-t2 
00049 10000000-1 
00059 31415926-1 
00069 18000000-3 
00079 20000000-1 
00089 90000000-2 
00099 50000000-1 
00109 85000000-2 
00119 30000000-1 
00129 14696000-2 
00139 OOOOOOOORR 
00143 0001 


00153 

T 1 

00163 

SSI 

00173 

G 

00183 

SHR 

002 13 

SIGD 

00223 

DELTAD 

002 33 

SIGR 

00243 

AX 

00273 

PP21 

00283 

UELTAR 

00293 

T2 

00303 

SS2 

00333 

THECR 

00343 

THECD 

00353 

THDN 

00363 

THDP 

00393 

PHIR 

00403 

TEST1 

00413 

TEST2 

00423 

TEST3 

00453 

WBL 

00463 

BFTWL 

00473 

CL 

00483 

DL 

00513 

FL 

00523 

FPL 

00533 

AWL 

00543 

8WL 



APPENDIX D 

SHEAR WAVE - SHOCK FRONT DATA 



M A = 

1 . 50 

DELTA = 

8.98 

SIGMA = 

54.43 

M = 1 . 16 THETA CR 

= 26.03 


THETA 

PH I —PR 

WS/WA 

DEL S 

WP/ WA 

DEL P 

(DP )RMS 

KSUR/K 

VSUR/A 

ADX/EL 

DEL SH 

-85 .0 

-81.98 

9.200 

0 .00 

-2.420 

0.00 

1.0163E+03 

.502 

2.408 

4.338 

90. 000 

-80.0 

-7 3.89 

4.39 7 

0.00 

-5.186 

0.00 

1 .087 1E+03 

.550 

2.301 

4.705 

90.000 

-75 .0 

-65.67 

2.643 

0.00 

-8.901 

0.00 

1.2406E+03 

.5 8,9 

2.230 

5.49 7 

90. 000 

-70.0 

-57.25 

1. 503 

0.00 

-15.251 

0.00 

1.5889E+03 

.618 

2.186 

7.2 84 

90.000 

-65.0 

-48 . 54 

-.022 

0 .00 

-32.881 

0.00 

2.7292E+03 

.637 

2.165 

' 13.100 

90.000 

-60.0 

-39.47 

33.790 

0.00 

518.270 

0.00 

3.5685E+04 

. 644 

2.169 

181.836 

90.000 

-55.0 

-29.92 

3.160 

0.00 

26. 140 

0.00 

1.5354E+03 

.637 

2.202 

8.457 

90.000 

-50.0 

-19.79 

2.092 

0.00 

11.030 

0.00 

8 . 643 8E +02 

.615 

2.273 

3.455 

90.000 

-45.0 

-8.92 

1. 607 

U .00 

4.707 

0.00 

2 . 1328E+02 

.575 

2.408 

1.55a 

90. 000 

-40.0 

-2.94 

1.215 

0.00 

-.549 

0.00 

2 .2369E+01 

.565 

2.405 

.094 

90.000 

-35 .0 

-16.22 

-1.064 

0.00 

-36. 638 

0.00 

1.3551E+03 

.744 

1.778 

10.137 

90. 000 

-30.0 

-32.13 

1.527 

0.00 

5.797 

0.00 

1 . 97 19E + 02 

.995 

1.284 

2.056 

90.000 

-25.0 

54.84 

1. 200 

-2.23 

1.134 

-41 .95 

3.5901E+01 

-.281 

-4.345 

.701 

69.082 

-20.0 

47.94 

1 .129 

-3.74 

1.197 

85.04 

3.5656E+01 

-.083 

-13*937 

.625 

42.908 

-15.0 

39.21 

1.085 

-3.44 

1.260 

56.87 

3.5 787E+0 1 

.115 

9.361 

.583 

28.899 

-10.0 

28.23 

1.057 

-2. 55 

1.312 

35.54 

3.5950E+01 

.313 

3.194 

.558 

18.148 

-5.0 

14.91 

1.041 

-1.35 

1. 344 

17.17 

3.6 063 E+0 1 

.509 

1.791 

.545 

8.792 

0.0 

0.00 

1 .036 

0.00 

1.355 

0.00 

3.6102E+01 

.701 

1.164 

. 540 

0.000 

5.0 

-14.91 

1.041 

1.35 

1.344 

-17.17 

3.6 063 E+0 1 - 

.887 

.804 

.545 

-8.792 

10.0 

-28.23 

1 .057 

2. 55 

1.312 

-35.54 

3 . 5950E+0 1 

1.067 

.567 

.558 

-18.148 

15 .0 

-39.21 

1.085 

3.44 

1.260 

-56.87 

3.5787E+01 

1.239 

.398 

.583 

-28.B99 

20.0 

-47 .94 

1 .129 

3.74 

1.197 

-85.04 

3.5656E+01 

1.401 

.268 

.625 

-42.908 

25.0 

-54.84 

1.200 

2.23 

1. 134 

41.95 

3 . 5 90 1 E+0 1 

1.552 

.165 

.701 

-69. 082 

30.0 

32.13 

1. 182 

0.00 

. 223 

0.00 

7.5914E-00 

.219 

.619 

.385 

-90.000 

35 .0 

16.22 

1. 205 

0.00 

.064 

0.00 

2.4038E-00 

.404 

.034 

.295 

-90. 000 

40.0 

2.94 

1.249 

0.00 

.003 

0.00 

1.3976E-01 

.5 09 

-.213 

.248 

-90.000 

45.0 

8.92 

1.314 

0.00 

-. 024 

0.00 

1 . 1258 E — 00 

.416 

-.552 

.218 

-90.000 

50.0 

19.79 

1.408 

0.00 

-.037 

0.00 

1.9 194E-O0 

.285 

-1.224 

. 196 

-90.000 

55 .0 

29.92 

1. 540 

0.00 

-.041 

0.00 

2 . 4509E-00 

.167 

-2.798 

.181 

-90. 000 

60.0 

3 9.47 

1.729 

0.00 

040 

0.00 

2.82 13E-00 

.057 

-10.141 

.169 

-90.000 

65.0 

48.54 

2.009 

0.00 

-.037 

0.00 

3 . 0849E-00 

-.044 

15.499 

.160 

-90. 000 

70.0 

57.25 

2.445 

0.00 

-.031 

0.00 

3.2732E-00 

-.139 

5.710 

.154 

-90.000 

75.0 

65.67 

3. 192 

0.00 

-.024 

0.00 

3 .405 IE-00 

-.226 

3.938 

. 149 

-90. 000 

BO .0 

73.89 

4.717 

0.00 

-.016 

0.00 

3 .492 5E-00 

-.306 

3.204 

.145 

-90.000 

85.0 

81.98 

9.349 

0 .00 

-.008 

0.00 

3 . 5423 E— 00 

-.379 

2.809 

. 144 

-90.000 


CD 



MA = 

2.00 

DELTA = 

8.96 

SIGMA = 

38.20 

M = 1.67 

THETA CR 

= 26. 

39 


theta 

P HI — PR 

WS/WA 

DEL S 

WP/WA 

DEL P 

{ DP )RMS 

KSUR/K 

VSUR/A 

ADX/EL 

DEL SH 

-85.0 

-82.09 

6.991 

0 .00 

-1.679 

0.00 

9.7343 E+0 2 

.382 

3.334 

4.032 

90.000 

-80.0 

-74.12 

3.336 

0.00 

-3.590 

0.00 

1 .0384E+03 

.449 

3.090 

4.362 

90.000 

-75.0 

-66.01 

1.999 

0 .oo 

-6.129 

0.00 

1 . 1782E+03 

.509 

2.927 

5.068 

90.000 

-70.0 

-57.68 

1.134 

0.00 

-10.369 

0.00 

1 .4889E+03 

.562 

2.816 

6.628 

90.000 

-65.0 

-49.07 

.025 

0 .00 

-21.435 

0.00 

2 . 4503 E+0 3 

.606 

2.741 

11.423 

90. 000 

-60.0 

-40.07 

-33.291 

0.00- 

•777.005 

0.00 

7.36 10E+04 

.641 

2.696 

364.472 

90.000 

-55.0 

-30 . 58 

2.612 

0 .00 

20.403 

0.00 

1 .6470E+03 

.666 

2.679 

8.820, 

90. 000 

-50.0 

-20.48 

1.663 

0.00 

8.322 

0.00 

5.8442E+02 

.678 

2.690 

3.480 

90.000 

-45.0 

-9.59 

1.260 

0.00 

3.599 

0.00 

2 . 2 346E+0 2 

.675 

2.740 

1.580 

90.000 

-40.0 

-2.33 

.949 

0.00 

-. 126 

0 .00 

7.0656E-00 

.683 

2.725 

.200 

90.000 

-35.0 

-15.79 

-.482 

0 .00 

-20.009 

0.00 

1 . 0 105 E+0 3 

.827 

2.247 

7.280 

90. 000 

-30. 0 

-32. 16 

1.187 

0.00 

4.092 

0.00 

1 .8970E+02 

1.021 

1.805 

1.946 

90.000 

-25.0 

53.39 

.925 

-2.80 

.855 • 

-48 .99 

3.6799E+01 

. 194 

9.314 

.682 

65.606 

-20.0 

46.41 

.870 

-4.15 

.9 04 

82.27 

3.6582E+01 

.337 

5.2 44 

. 608 

41.594 

-15.0 

37.7 2 

.835 

-3.78 

.955 

55.18 

3.6744E+01 

.477 

3.585 

.566 

28. 119 

-L0.0 

26.97 

.813 

-2.79 

.995 

34.53 

3.6929E+01 

.614 

2.673 

.542 

17.686 

-5.0 

14. 17 

.801 

-1.47 

1.021 

16.69 

3. 7 0 54 E +01 

.746 

2.089 

.529 

8.575 

0.0 

0.00 

.796 

0.00 

1.030 

0.00 

3 . 7097E+01 

.872 

1.678 

.525 

0.000 

5.0 

-14. 17 

.801 

1.47 

1.021 ■ 

-16.69 

3 . 7054E+0 1 

.992 

1.369 

.52 9 

-8.575 

10.0 

-26.97 

.813 

2.79 

.993 ■ 

-34.53 

3.6929E+01 

1.104 

1.125 

.542 

-17.686 

15.0 

-37 .72 

.835 

3.78 

.955 -55.18 

3.6744E+0 1 

1.207 

.924 

. 566 

-28. 119 

20.0 

-46.41 

.870 

4.15 

.904 • 

-82.27 

3.6582E+01 

1.302 

.754 

.608 

-41.594 

25.0 

-53. 39 

.925 

2.80 

.855 

48.99 

3 . 6 799E+0 1 

1.386 

.605 

.682 

-65.606 

30.0 

32. 16 

.912 

0.00 

.179 

0.00 

8.3 17 9E-00 

.48 9 

1.413 

.389 

-90.000 

35.0 

15.79 

.926 

0.00 

.049 

0.00 

2.5 104E-00 

.601 

.895 

.296 

-90.000 

40.0 

2. 33 

.958 

0.00 

0.000 

0.00 

5.0358E-02 

.653 

.583 

.247 

-90.000 

45.0 

9.59 

1.006 

0.00 

-.021 

0.00 

1.3103 E -00 

.558 

.394 

.217 

-90. 000 

30.0 

20.48 

1.076 

0.00 

-.030 

0.00 

2 . 158 9E -00 

.443 

. 131 

. 195 

-90.000 

55.0 

30.58 

1.175 

0 .00 

-.033 

0.00 

2 . 72 54 E— 00 

.334 

-.311 

. 180 

-90. 000 

60.0 

40.07 

1 . 318 

0.00 

-.032 

0.00 

3. 1193E-00 

.230 

-1.152 

.168 

-90.000 

65.0 

49.07 

1.530 

0.00 

029 

0 .00 

3 .3 993 E— 00 

. 130 

-3.258 

.159 

-90. 000 

70.0 

57 .68 

1.861 

0.00 

-.025 

0.00 

3 . 5990E-00 

.034 

-16.995 

.152 

-90.000 

75 .0 

66.01 

2.428 

0.00 

-.019 

0.00 

3.7389E-00 

-.058 

12.601 

. 148 

-90. 000 

80.0 

74. 12 

3.58 7 

0.00 

-.013 

0.00 

3.83 1 5E-00 

-.146 

6.005 

.. 144 

-90.000 

85.0 

82.09 

7.108 

0 .00 

-. 006 

0.00 

3.8843E-00 

-.230 

4.427 

. 143 

-90. 000 



MA » 

2.50 

DELTA = 

8.96 

SIGMA = 

30.88 

M = 2.12 

THETA CR 

= 27. 

20 


THETA 

PHI— PR 

WS/WA 

DEL S 

WP/WA 

DEL P 

{DP )RMS 

KSUR/K 

VS UR / A 

ADX/EL 

DEL SH 

-85.0 

-82. 39 

5.795 

0.00 

-1. 112 

0.00 

8.8393E+02 

.324 

4. 157 

3.376 

90.000 

-80.0 

-74.71 

2.754 

0.00 

-2.363 

0.00 

9.3677E+02 

.398 

3.782 

3.629 

90.000 

-75.0 

-66 .88 

1. 640 

0.00 

-3.982 

0.00 

1 .0480E+03 

.466 

3.541 

4.159 

y 0. 000 

-70.0 

-58.82 

.925 

0.00 

-6. 542 

0.00 

1 .2842E+03 

.528 

3.379 

5.277 

90.000 

-65 .0 

-50.45 

.096 

0.00 

-12.369 

0.00 

1 .9291E+03 

.583 

3.266 

8.307 

90.000 

-60.0 

-41.67 

-4. 574 

0.00 

-59,411 

0.00 

7 . 6605E + 03 

.629 

3. 191 

35.068 

90.000 

-55.0 

-32.35 

2.832 

0.00 

19.287 

0.00 

2. 1128E+03 

.667 

3. 145 

10.473 

90. 000 

-50.0 

-22.34 

1.570 

0.00 

6.928 

0.00 

6.5796E+02 

.694 

3. 128 

3.630 

90.000 

-45.0 

-11.45 

1.143 

0.00 

3.048 

0.00 

2.5497E+02 

.709 

3. 145 

1.662 

90. 000 

-40.0 

-.64 

.850 

0.00 

.389 

0.00 

2.9058E+01 

.713 

3. 182 

.427 

90.000 

-35.0 

-14.52 

.020 

0.00 

-8.022 

0.00 

5 .4055E+02 

.839 

2.734 

3.478 

90. 000 

-30.0 

-32. 10 

1.045 

0.00 

2.854 

0.00 

1.7557E+02 

1.006 

2.286 

1.717 

90.000 

-25 .0 

49.8 8 

.796 

-4.20 

.692 

-62.00 

3 . 9335E+0 1 

.439 

5.211 

.640 

59.259 

-20.0 

42.81 

. 746 

-5.27 

. 738 

76.26 

3 .9208E+01 

.546 

4. 132 

.569 

38.800 

-15 .0 

34.29 

.715 

-4.70 

.785 

51 .44 

3.94666+01 

.650 

3. 403 

.530 

26.427 

-10.0 

24. 17 

.69 5 

-3.46 

.823 

32.26 

3.9717E+01 

.748 

2.871 

. 508 

16.677 

-5 .0 

12. 55 

.684 

-1.82 

.847 

15.61 

3. 9880E+01 

.841 

2.460 

.495 

8.098 

0.0 

0.00 

.681 

0.00 

.855 

0.00 

3 .9935E+01 

.927 

2. 128 

.492 

0.000 

5.0 

-12.55 

.684 

1.82 

.847 

-15.61 

3.9880E+01 

1.007 

1.851 

.495 

-8.098 

10.0 

-24. 17 

.695 

3.46 

.823 

-32.26 

3.9717E+01 

2.078 

1. 612 

.508 

-16.677 

15.0 

-34.29 

.715 

4.70 

.785 

-51 .44 

3 • 9466E+0 1 

1.142 

1.40 2 

. 530 

-26.427 

20.0 

-42.81 

. 746 

5.27 

.738 

-76.26 

3 . 9208E+0 1 

1 . 196 

1.212 

.569 

-38.800 

25.0 

-49.88 

. 796 

4.20 

.692 

62.00 

3.9335E+01 

1.242 

1.038 

. 640 

-59.259 

30.0 

32.10 

.788 

0.00 

. 168 

0.00 

1.0383E+01 

.600 

1.864 

.401 

-90.000 

35.0 

14. 52 

.791 

0.00 

.040 

0.00 

2.7 172E-00 

.680 

1.38 1 

.296 

-90. 000 

40.0 

.64 

.813 

0.00 

-.003 

0.00 

2.8176E-01 

.707 

1.064 

.245 

-90.000 

45.0 

11.45 

.849 

0.00 

-.022 

0.00 

1 .8 980E-00 

.602 

.931 

.214 

-90.000 

50.0 

22.34 

.904 

0.00 

-.030 

0.00 

2.8929E-00 

.497 

.732 

. 192 

-90.000 

55 .0 

32.35 

.985 

0 .00 

-.032 

0.00 

3.5520E-00 

.396 

.416 

.177 

-90. 000 

60.0 

41.67 

1.102 

0.00 

-.031 

0.00 

4.008 IE-00 

- .297 

119 

.165 

-90.000 

65 .0 

50.45 

1.276 

0.00 

-.027 

0.00 

4.33 12E-00 

.200 

-1. 174 

.156 

-90. 000 

70.0 

58.82 

1.549 

0.00 

-.023 

0.00 

4.5610E-00 

.106 

-4.090 

.150 

-90.000 

75.0 

66.88 

2.019 

0.00 

-.017 

0.00 

4.72 18 E— 00 

.013 

-45.677 

.145 

-90. 000 

80.0 

74.71 

2.979 

0.00 

-.012 

0.00 

4.828 IE-00 

-.076 

10.790 

. 142 

-90.000 

85.0 

82. 39 

5.901 

0.00 

-.006 

0.00 

4.8887E-00 

162 

6. 174 

.140 

-90.000 


CO 

CO 



MA = 

3.00 

DELTA = 

8.97 

SIGMA = 

26.46 

M » 2.55 THETA CR 

= 27.83 



THETA 

( 

PHI-PR 

WS/WA 

DEL S 

WP/WA 

DEL P 

(DP )RMS 

KSUR/K 

VSUR/A 

ADX/EL 

DEL SH 

-85.0 

-82.70 

5.026 

0.00 

-.773 

0.00 

8. 1665E+02 

.287 

4.940 

2.850 

90. 000 

-80.0 

-75. 32 

2.381 

0.00 

-1.634 

0.00 

8.5999E+02 

. 364 

4.433 

3.044 

90.000 

-75.0 

-67.80 

1.409 

0 .00 

-2.720 

0.00 

9.4958E+02 

.437 

4. 121 

3.444 

90.000 

-70.0 

-60.03 

.794 

0.00 

-4.353 

0.00 

1 . 1317E+03 

. 504 

3.912 

4.252 

90.000 

-65.0 

-51.92 

.136 

0.00 

-7.666 

0.00 

1 . 58 06E+0 3 

.565 

3.767 

6.227 

90. 000 

-60.0 

-43.38 

-1.916 

0.00 

-22.553 

0.00 

3.8352E+03 

.618 

3.666 

16.076 

90 . 000 

-55.0 

-34.25 

3.707 

0.00 

22.214 

0.00 

3. 1999E+03 

. 66 4 

3.599 

14.537 

90. 000 

-50.0 

-24. 38 

1.600 

0.00 

6.263 

0.00 

7.7940E+02 

.700 

3.562 

3.943 

90.000 

-45.0 

-13. 53 

1. 102 

0.00 

2.732 

0.00 

2.9821E+02 

.725 

3.557 

1.773 

90. 000 

-40.0 

-1. 32 

.807 

0.00 

.663 

0.00 

6.4323E+01 

.735 

3.596 

.615 

90.000 

-35.0, 

-12.94 

.245 

0.00 

-3.582 

0.00 

3. 1180E+02 

.8 37 

3.213 

1.718 

90. 000 

-30.0 

-31.80 

.966 

0.00 

2.130 

0.00 

1.6831E+02 

.987 

2.751 

1.554 

90.000 

-25.0 

46.62 

.717 

-5.65 

.584 

-71.24 

4.2 34 1E+0 1 

.576 

4.723 

.60 5 

54. 859 

-20.0 

39. 55 

.669 

-6. 52 

.628 

71.32 

4 . 2340E+0 1 

.663 

4.079 

.536 

36.602 

-15.0 

31.30 

.640 

-5.76 

.‘674 

48.28 

4 . 2 724E+0 1 

.744 

3.582 

.499 

25.068 

-10.0 

21.81 

.622 

-4.22 

.710 

30.33 

4.3062E+01 

.820 

3. 180 

.478 

15.859 

-5.0 

11.23 

.612 

-2.22 

.734 

14.68 

4.3275E+01 

.890 

2.8 44 

.466 

7.710 

0.0 

0.00 

.609 

0.00 

.742 

0.00 

4 . 3346E+0 1 

.953 

2.555 

.462 

0.000 

5.0 

-11.23 

.612 

2.22 

.734 

-14.68 

4.3275E+01 

1.009 

2. 299 

.466 

-7. 710 

10.0 

-21.81 

.622 

4. 22 

.710 

-30.33 

4.3062E+01 

1.057 

2.069 

.478 

-15.859 

15.0 

-31.30 

. 640 

5.76 

.674 

-48.28 

4.2724E+01 

1.097 

1.858 

.49 9 

-25.068 

20.0 

-39,55 

.669 

6. 52 

.628 

-71.32 

4.2340E+01 

1.129 

1.660 

.536 

-36.602 

25.0 

- 46 . 6 2 

.717 

5.65 

.584 

71.24 

4 . 2 341E+0 1 

1. 152 

1.471 

.60 5 

-54.859 

30.0 

31.80 

.714 

0.00 

.161 

0.00 

1 . 2769E+0 1 

« 664 

2.262 

.412 

-90.000 

35.0 

12.94 

.707 

0.00 

.031 

0.00 

2 . 7 772E-00 

.724 

1.794 

.29 5 

-90. 000 

40.0 

1.32 

.721 

0.00 

-.008 

0.00 

8.2032E-01 

.725 

1.498 

.242 

-90. 000 

45.0 

13.53 

.750 

0.00 

-.024 

0.00 

2.7133E-00 

.623 

1.388 

.210 

-90. 000 

50.0 

24. 38 

.795 

0.00 

-.031 

0.00 

3.865 IE-00 

.525 

1.212 

. 189 

-90.000 

55.0 

34.25 

.863 

0.00 

-.032 

0.00 

4 . 623 IE-00 

.429 

.940 

.173 

-90. 000 

60.0 

43.38 

.962 

0.00 

-.030 

0.00 

5. 1453E-00 

. 334 

.500 

.161 

-90.000 

65.0 

51.92 

1.112 

0.00 

-. 026 

0.00 

5.5 141E-00 

.240 

-. 289 

.153 

-90. 000 

70.0 

60.03 

1.348 

0.00 

-.022 

0.00 

5.77606-00 

.147 

-2.072 

. 146 

-90.000 

75.0 

67.80 

1. 755 

0.00 

-.017 

0.00 

5.9589E-00 

.056 

-9.617 

. 142 

-90. 000 

80.0 

75. 32 

2.588 

0.00 

-.011 

0.00 

6.07976-00 

-.033 

22.908 

.138 

-90.000 

85.0 

82.70 

5.124 

0.00 

-.005 

0.00 

6. 1486E-00 

-.121 

8.207 

.137 

-90.000 


100 



MA = 

1.75- 

DELTA = 

16. 73 

SIGMA = 

57.50 

M = 1.08 

THETA CR 

= 28. 

59 


THETA 

PHI -PR 

WS/WA 

DEL S 

WP/ WA 

DEL P 

(DP )RKS 

ksur/k 

VSUR/A 

AUX/EL 

DEL SH 

—85 .0 

-83.37 

9.494 

0.00 

-.935 

0.00 

7.3278E+02 

.555 

2.44 5 

2.056 

90. 000 

-80. 0 

-76 .66 

4.474 

0.00 

-1.956 

0.00 

7.6269E+02 

.609 

2.321 

2.171 

90.000 

-75.0 

-69.79 

2.630 

0.00 

—3.1 85 

0.00 

8.2223E+02 

.655 

2.230 

2.400 

90. 000 

-70.0 

-62.66 

1.490 

0.00 

-4.875 

0.00 

9 .346 1E+02 

.691 

2.164 

2.827 

90.000 

-65.0 

-55. 17 

.418 

0.00 

-7.724 

0.00 

1 . 16 99E+03 

.716 

2.119 

3.713 

90. 000 

-60.0 

-47.20 

-1.456 

0.00 

-14.924 

0.00 

1 .8553E+03 

.731 

2.0 93 

6.271 

90.000 

-55 .0 

-38.58 

-32.856 

0.00- 

•155.549 

0.00 

1.6284E+04 

.733 

2.08 7 

59.705 

90.000 

-50.0 

-29.12 

5.017 

0.00 

1 5. 592 

0.00 

1 . 39996+03 

.720 

2.108 

5.714 

90.000 

-45.0 

-18.51 

2. 777 

0.00 

5.852 

o.'oo 

4.5680E+02 

.690 

2.167 

2. 171 

90.000 

-40.0 

-6.29 

1.915 

0.00 

2.171 

0.00 

1.48986+02 

.635 

2.299 

.945 

90.000 

-35.0 

-8.51 

1.082 

0 .00 

-1.483 

0.00 

9 . 03 53E+0 1 

.700 

2.021 

. 197 

90i 000 

-30.0 

-29.74 

2.166 

0.00 

3.404 

0.00 

1.8558E+02 

.979 

1.388 

1. 436 

90.000 

-25 .0 

40.76 

1.467 

-8.85 

1.023 

-83.83 

5 . 0448E+0 1 

. 176 

7.329 

.547 

49.261 

-20.0 

33.93 

1.357 

-9.52 

1 . 129 

63.65 

5 .09016 + 01 

.298 

4.067 

.482 

33.534 

-15.0 

26.3 5 

1. 292 

-8.35 

1.236 

43.18 

5 . 1709E+01 

.419 

2.695 

.448 

23. 131 

-10.0 

18.05 

.1.252 

-6.15 

1 . 325 

27.15 

5.2328E+01 

.536 

1.930 

.428 

14.682 

-5.0 

9.18 

1. 230 

-3.25 

1.384 

13.15 

5.2698E+01 

.649 

1.436 

.417 

7. 149 

0.0 

0.00 

1.223 

0.00 

1.405 

0.00 

5.2819E+01 

.757 

1.087 

.414 

0.000 

5 .0 

-9.18 

1.2 30 

3.25 

1.384 

-13.15 

5.26986+01 

.859 

.823 

.417 

-7. 149 

10.0 

-18.05 

1.2 52 

6.15 

1.325 

-27.15 

5 . 2328E+0 1 

.955 

.613 

.428 

-14.682 

15 .0 

-26.35 

1. 292 

8.3 5 

1.236 

-43.18 

5. 1709E+0 1 

1.044 

.441 

. 448 

-23. 131 

20.0 

-33.93 

1.357 

9.52 

1.129 

-63 .65 

5 . 090 IE +0 1 

1 . 124 

.294 

.482 

-33.534 

25.0 

-40 .76 

1.467 

8.85 

1.023 

83.83 

5 . 0448E+0 1 

1.196 

. 167 

.547 

-49.261 

30.0 

29.74 

1.483 

0.00 

. 327 

0.00 

1.78646+01 

. 332 

.200 

.425 

-90.000 

35.0 

8.51 

1.419 

0.00 

.032 

0.00 

1 .9686E-00 

.540 

-.123 

.2 87 

-90.000 

40.0 

6.29 

1.423 

0.00 

-.043 

0.00 

2.96966-00 

.524 

-.380 

.233 

-90.000 

45.0 

18.51 

1.462 

0.00 

070 

0.00 

5 . 4759E-00 

.380 

-.870 

.201 

-90. 000 

50.0 

29.12 

1.538 

0.00 

-.077 

0.00 

6 . 975 2E-00 

.253 

-1.821 

. 180 

-90.000 

55.0 

38.58 

1.658 

0 .00 

-.075 

0.00 

7 « 9522E-00 

. 135 

-4. 324 

.165 

-90. 000 

60.0 

47.20 

1.841 

0.00 

-.069 

0.00 

8.62126-00 

.026 

-27.102 

.153 

-90.000 

65 .0 

55. 17 

2. 119 

0.00 

-.060 

0.00 

9.09166-00 

-.076 

10.749 

. 145 

-90. 000 

70.0 

62.66 

2.562 

0.00 

-.049 

0.00 

9 .4248E-00 

-. 172 

5.393 

. 139 

-90.000 

75 .0 

69.79 

3.329 

0.00 

-.037 

0.00 

9.6569E-00 

-.263 

3.935 

. 134 

-90. 000 

60.0 

76.66 

4.903 

0.00 

--.025 

0.00 

9.8 10 IE-00 

-.346 

3.257 

.131 

-90.000 

85 .0 

83.37 

9.699 

0.00 

-.012 

0.00 

9.8973E-00 

-.423 

2.869 

. 129 

-90. 000 


TOT 



MA = 

2.00 

DELTA = 

16. 73 

SIGMA = 

47 .80 

M = 1.37 ' 

THETA CR 

= 28.60 



THETA 

PHI -PR 

WS/WA 

DEL S 

WP/WA 

DEL P 

(DP )RMS 

KSUR/K 

VSUR/A 

ADX/EL 

DEL SH 

-85 .0 

-83.39 

8.339 

0.00 

-.809 

0.00 

7 .3098E+02 

.463 

3.007 

2.034 

90.000 

-80.0 

-76.71 

3.929 

0.00 

-1.691 

0 .00 

7 .6056E+02 

.528 

2.801 

2.147 

90.000 

-75.0 

-69.86 

2. 309 

0.00 

-2.752 

0.00 

8. 1932E+02 

.585 

2.654 

2.371 

90. 000 

-70.0 

-62.75 

1 . 309 

0.00 

-4. 206 

0.00 

9 . 299 6E +02 

.635 

2.546 

2.789 

90.000 

-65.0 

-55.29 

.371 

0 .00 

— 6. 645 

0.00 

1 . 1605E+03 

.676 

2.467 

3.653 

90.000 

-60.0 

-47 . 3 A 

-1.245 

0.00 

-12.725 

0.00 

1 .8238E+03 

.708 

2.413 

6.113 

90.000 

-55.0 

-38 .74 

-23.634 

0 .00- 

■111.049 

0.00 

1.3399E+04 

.728 

2. 38 0 

48.718 

90.000 

-50.0 

-29.30 

4.517 

,0.00 

13.932 

0.00 

1 .4413E+03 

.736 

2.371 

5.834 

90.000 

-45.0 

-18.70 

2. 469 

0 .00 

5. 172 

0.00 

4.6515E+02 

.729 

2.39 5 

2.192 

90. 000 

-40.0 

-6.49 

1.697 

o.Oo 

1.931 

0.00 

1 .5259E+02 

.701 

2.470 

.957 

90.000 

-35 .0 

-8.32 

.968 

0 .00 

-1.216 

0.00 

8.5279E+01 

.763 

2.233 

. 167 

90. 000 

-30.0 

-29.62 

1.918 

0.00 

2.989 

0.00 

1 .8752E+02 

.995 

1.671 

1.439 

90.000 

-25.0 

40.57 

1. 292 

-8.97 

.896 

-84.16 

5.0788E+01 

.375 

4. 290 

.546 

49. 128 

-20.0 

33.76 

1 . 194 

-9.64 

. 989 

63 .43 

5. 12 64E+01 

.48 0 

3.217 

.480 

33.458 

-15 .0 

26.20 

1.138 

-8.45 

1.084 

43.03 

5.2093E+01 

.582 

2.525 

. 446 

23-082 

-10.0 

17.94 

1 . 102 

-6.23 

1.163 

27.05 

5.2726E+01 

.678 

2.035 

.426 

14.652 

-5 .0 

9.12 

1.083 

-3.30 

1.216 

13.10 

5.3102E+01 

.770 

1 • 664 

.416 

7. 134 

0.0 

0.00 

1.076 

0.00 

1.234 

0.00 

5 . 3226E+01 

.856 

1.371 

.413 

0.000 

5.0 

-9.12 

1.083 

3.30 

1.216 

-13.10 

5.3102E+01 

.935 

1. 129 

.416 

-7. 134 

10.0 

-17.94 

1 . 102 

6.23 

1.163 

-27.05 

5.2726E+01 

1.008 

.924 

.426 

-14.652 

15 .0 

-26 .20 

1.138 

8.45 

1.084 

-43.03 

5.2093E+01 

1.072 

.745 

.446 

-23. 082 

20.0 

-33.76 

1 . 194 

9.64 

.989 

-63.43 

5. 12 64E+01 

1.129 

.585 

.480 

-33.458 

25.0 

-40.57 

1.292 

8.97 

.896 

84.16 

5.0788E+0L 

1.176 

.439 

.546 

-49. 128 

30.0 

29.62 

1.307 

0.00 

.287 

0.00 

1.8013E+01 

.487 

.757 

.426 

-90.000 

35.0 

8.32 

1. 249 

0.00 

.027 

0.00 

1 .9059E-00 

.639 

.342 

.287 

-90.000 

40.0 

6 . 49 

1.251 

0.00 

-.038 

0.00 

3 . 07 98 E-00 

.611 

. 109 

.233 

-90.000 

45.0 

18.70 

1.286 

0.00 

-. 062 

0.00 

5 .608 IE-00 

.48 2 

-. 177 

.201 

-90.000 

50.0 

29.30 

1.352 

0.00 

-.068 

0.00 

7 . 1199E-00 

. 364 

-.651 

.179 

-90.000 

55.0 

38.74 

1.457 

0 .00 

-.067 

0.00 

8. 1049E-00 

.253 

-1.526 

. 164 

-90.000 

60.0 

47.34 

1.617 

0.00 

-.061 

0.00 

8 . 7792E-00 

. 148 

-3.606 

.153 

-90.000 

65.0 

55.29 

1.862 

0.00 

-.052 

0.00 

9.2534E-00 

.047 

-14. 386 

. 145 

-90. 000 

70.0 

62.75 

2.2 50 

0.00 

-.043 

0.00 

9.589 IE — 00 

-.049 

16.355 

.138 

-90.000 

75 .0 

69.86 

2.924 

0.00 

-.033 

0.00 

9.8231E-O0 

-. 142 

6 • 644 

. 134 

-90. 000 

80.0 

76.71 

4.306 

0.00 

-.022 

0.00 

9.977 4E-00 

-.230 

4.646 

.131 

-90.000 

85 .0 

83. 39 

8.519 

0 .00 

-.on 

0.00 

1 .0065E+01 

-.313 

3.786 

. 129 

-90. 000 


102 




38.90 

M = 1.79 THETA CR 

= 28.70 


DEL P 

(DP)RMS 

KSUR/K 

VSUR/A 

ADX/EL 

06L SH 

0.00 

7. 1355E+02 

.380 

3.897 

1.747 

90. 000 

0.00 

7 . 3861E+02 

.453 

3.550 

1.835 

90.000 

0.00 

7.8759E+02 

.521 

3.312 

2.006 

90. 000 

0.00 

8.7702E+02 

.582 

3.141 

2.315 

90.000 

0.00 

1.0522E+03 

.636 

3.014 

2.915 

90.000 

0.00 

1.4857E+03 

.682 

2.920 

4.384 

90.000 

0.00 

3.88076+03 

.719 

2.853 

12.420 

90.000 

0.00 

2.47 146+03 

.746 

2.811 

8.798 

90.000 

0.00 

6.22 14E+02 

.760 

2.796 

2.567 

90. 000 

0.00 

2.132 2E +02 

.757 

2.821 

1.130 

90.000 

0.00 

2.06596+01 

.787 

2.706 

. 181 

90.000 

0.00 

2.3728E+02 

.971 

2.174 

1.581 

90.000 

-88.20 

5.6182E+01 

.571 

3.633 

.524 

47. 614 

60.59 

5.7065E+01 

.653 

3.096 

.459 

32.570 

41.07 

5.8248E+01 

.730 

2.67 6 

.426 

22.510 

25.80 

5.9108E+01 

.801 

2.334 

.407 

14.300 

12.49 

5.9609E+01 

.867 

2.045 

. 397 

6.965 

0.00 

5.9772E+01 

.926 

1.796 

. 393 

0.000 

-12.49 

5 .9609E+0 1 

.977 

1.57 5 

.397 

—6 .965 

-25.80 

5.9108E+01 

1.022 

1.374 

.407 

-14.300 

-41.07 

5 . 8 248E+0 1 

1.058 

1. 189 

. 426 

-22. 510 

-60.59 

5 . 7065E+0 1 

1.0,8 7 

1.015 

.459 

-32.570 

88.20 

5.6182E+01 

1.107 

.849 

.524 

-47. 614 

0.00 

1 .95456+01 

.632 

1.215 

.422 

-90.000 

0.00 

6 . 3790E-0 1 

.729 

.812 

.280 

-90.000 

0.00 

5.0137E-00 

.661 

.622 

.2 26' 

-90.000 

0.00 

7 .B562E-00 

.549 

.413 

. 195 

-90. 000 

0.00 

9 . 5494E-00 

.444 

.091 

.174 

-90.000 

0.00 

1.0650E+01 

.342 

-.425 

.159 

-90.000 

0.00 

1. 1402E+01 

.243 

-1.360 

. 148 

-90.000 

0.00 

1.1931E+01 

.145 

-3.520 

. 140 

-90.000 

0.00 

1 .23O5E+01 

.050 

-13.657 

. 134 

-90.000 

0.00 

1 .2565E+01 

-.042 

20.565 

.130 

-90. 000 

0.00 

1 .27376+01 

-.132 

7.817 

. 127 

-90.000 

0.00 

1 .2835E+01 

-.218 

5.446 

. 125 

-90.000 


103 



MA = 

3.00 

DELTA = 

16.75 

SIGMA = 

34.09 

M = 2.16 THETA CR 

= 28.63 



THETA 

P HI —PR 

WS/WA 

DEL S 

WP/WA 

DEL P 

(DP ) RMS 

KSUR/K 

VSUR/A 

ADX/EL 

DEL SH 

-85.0 

-84. 10 

6.303 

0.00 

-.390 

0.00 

7. 1117E+02 

.334 

4.699 

1.497 

90.000 

-80.0 

-78.13 

2.961 

0.00 

-.8 09 

0.00 

7.32 53E+02 

.412 

4.214 

1.5 64 

90.000 

-75.0 

-71.99 

1.737 

0.00 

-1.293 

0.00 

7.7362E+02 

.48 4 

3.893 

1.694 

90. 000 

-70.0 

-65.60 

.999 

0.00 

-1.909 

0.00 

8 .465 0E +02 

.551 

3.666 

1.921 

90.000 

-65 .0 

-58.83 

.369 

0.00 

-2.812 

0.00 

9.8 173E+02 

.611 

3.500 

2. 338 

, 90.000 

-60.0 

-51.58 

-.440 

0.00 

-4.482 

0.00 

1.2783E+03 

.665 

3.374 

3.242 

90.000 

-55 .0 

-43.65 

-2. 530 

0.00 

-9.712 

0.00 

2.3 191E+03 

.710 

3.28 0 

6.376 

90. 000 

-50.0 

-34.82 

16.606 

0.00 

41.191 

0.00 

8.3749E+03 

.7 46 

3.213 

25.581 

90.000 

-45.0 

-24.75 

3.007 

0.00 

5.271 

0.00 

9.2363E+02 

.77 2 

3. 173 

3.256 

90. 000 

-40.0 

-12.90 

1.756 

0.00 

2.010 

0.00 

3.0633E+02 

.783 

3.167 

1.353 

90.000 

-35.0 

-1.86 

1. 114 

0.00 

.314 

0.00 

4. 1903E+01 

.789 

3. 158 

.447 

90.000 

-30.0 

-23.93 

2.377 

0.00 

3.620 

0.00 

4.2565E+02 

.941 

2.644 

2.294 

90.000 

-25 .0 

35.80 

1.060 

-12.44 

.605 

88.70 

6.31896+01 

.670 

3. 676 

. 506 

46.744 

-20.0 

29.38 

.968 

-13.25 

.692 

58.13 

6.4733E+01 

.739 

3.274 

.440 

32.033 

-15.0 

22.51 

.916 

-11.78 

.783 

39.30 

6 • 6464E+0 1 

.802 

2.937 

.40 7 

22. 156 

-10.0 

15.25 

.885 

— 8.81 

.861 

24.66 

6.7666E+01 

.859 

2.648 

.389 

14.078 

-5 .0 

7.70 

.868 

-4.72 

.914 

11.93 

6.8351E+01 

.910 

2.394 

.379 

6.857 

0.0 

0.00 

.862 

0.00 

.934 

0.00 

6.8573E+01 

.954 
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